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MANUSCRITS RECUS LE 31 JANVIER 1961 
Première Section 


SCIENCES MATHEMATIQUES ET ASTRONOMIQUES 


Systèmes matériels de trois masses ponctuelles 
MOUVEMENT PLAN EN COORDONNÉES SYMÉTRIQUES 


PAR 


ANDRÉ DEPRIT et ANDRÉ DELIE 


Faculté des Sciences 
Université de Louvain 


Dans le cas particulier du problème plan des trois corps, F.D. 
Murnaghan a donné la fonction hamiltonienne et les équations cano- 
niques du mouvement au moyen de quatre coordonnées symétriques, 
parmi lesquelles les distances mutuelles entre les trois points ($ 1). 

Cette fonction hamiltonienne contient une forme quadratique 
non diagonalisée dans les moments; en outre les équations canoniques 
qu’on en dérive ne sont pas régularisées pour les chocs doubles. 

Mettant a profit des coordonnées périmétriques introduites par 
Monseigneur Lemaitre, nous transformons la fonction hamiltonienne 
de F.D. Murnaghan et nous obtenons une fonction régularisée pour les 
chocs binaires (§ 2). 

Ensuite, par des transformations élémentaires (passage a des 
coordonnées polaires sphériques et projection stéréographique), nous 
arrivons a une fonction lagrangienne et une fonction hamiltonienne 
associée qui nous semble marquer des avantages décisifs sur les fonctions 
correspondantes de F.D. Murnaghan : 

1°) la forme quadratique dans les vitesses et la forme adjointe 
correspondante sont diagonalisées; 

2°) la fonction hamiltonienne est régularisée pour les chocs binaires; 

30) le moment polaire figure parmi les variables canoniques, ce 
qui équivaut à faire figurer l'identité de Lagrange parmi les équations 
du mouvement, une fois la réduction faite par l’intégrale des aires; 

4°) la fonction hamiltonienne s’adresse à toute fonction de forces 


qui ne dépend que des distances, et pas seulement à la seule fonction 
des forces gravifiques. (§3 et 4). 

Les calculs littéraux sont fort longs; si nous avions dû, dans chaque 
cas, indiquer la meilleure maniére d’aboutir aux expressions sous la 
forme la plus symétrique et la plus simplifiée que nous avons cru pou- 
voir leur donner, nous aurions augmenté démesurément le volume de 
cet article. 

A propos des réductions que nous effectuons ici pour le probleme 
des trois corps, nous avons eu de nombreuses conversations avec Mon- 
seigneur Lemaitre; nous l’en remerçions vivement. 

Mademoiselle Hélène De Laet a bien voulu vérifier avec nous 
quelques-uns de nos calculs les plus lourds. 

Ce travail reprend à sa manière l’équivalence établie par Mon- 
sieur l'Abbé Steenbeckeliers et Madame Deprit-Bartholomé entre la 
fonction hamiltonienne en coordonnées symétriques de A. Wintner 
et les fonctions hamiltoniennes de Monseigneur Lemaitre. 

Soient trois points P; affectés chacun respectivement d’une masse 
mi > 0. On suppose que les trois points demeurent dans un plan fixe 
et qu'ils sont soumis à des forces dérivables d’une fonction de forces U 
qui ne dépend que des distances mutuelles entre les points. 


C’est pourquoi le moment cinétique € du système est un vecteur 
constant perpendiculaire au plan du mouvement et le centre de gravité G 
se déplace d’un mouvement rectiligne uniforme. 

Dans le plan du mouvement, on choisit un repère barycentrique 
Gxy de telle manière que l'orientation définie dans le plan par ce repère 
induise sur le triangle des masses l'orientation P;P2P3. 

Pour des raisons de concision, nous conviendrons que, toutes les 
fois que, dans une relation, une définition ou un énoncé, le triplet 
(i,j,k) figure au complet, il sera entendu, sauf mention explicite du 
contraire, qu’il parcourt l’ensemble des permutations paires effectuées 
sur l’ensemble (1, 2, 3). 
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Deuxième Section 


SCIENCES PHYSIQUES ET CHIMIQUES 


Sur les aberrations périaxiales de l'œil. I 


PAR 
A. BIOT 


SOMMAIRE. Supposant que l’indice du cristallin, dans l’œil théorique de 
Y. Le Grand, varie avec la distance à l’axe de figure, on examine une méthode 
de calcul simplifiée des aberrations et diverses lois de variation. On est amené 
à étudier une loi de la forme simple n° = n + by? dans laquelle n, est l’indice 
sur l’axe. On établit l’équation numérique du rayon courbe dans le cristallin 
pour 6 hauteurs d’incidence sur la cornée de rayons parallèles faisant un 
angle d’environ 5° avec l’axe. On calcule ensuite les hauteurs d’intersection 
de ces rayons avec le plan focal image. Le minimum de dispersion de ces 
points a lieu pour b voisin de + 0,000668. Il y a alors une correction notable 
des mêmes aberrations calculées pour l’œil théorique ordinaire. Cette cor- 
rection s'étend même au-delà de la zone de diamètre 4 mm. 


1. Le problème. Dans un travail précédent (1) nous avons précisé 
et calculé les aberrations géométriques que l’on rencontre en la fovéa 
de l’œil théorique de Y. Le Grand (non accommodé). Nous avons en 
particulier recherché les hauteurs d’intersection, dans le plan focal 
image, de rayons incidents qui devraient se rencontrer en la fovéa 
s’il n’y avait pas d’aberration. Celles que l’on trouve ainsi sont plus 
grandes que celles que l’on observe sur l’œil réel. Supposant que l’in- 
dice, dans le cristallin, varie avec la distance à l’axe de figure, nous 
recherchons une méthode de calcul et une loi de variation compatible 


avec les observations. 


2. Etude préliminaire. a) Le même problème a été étudié en sup- 
posant que les rayons lumineux restent rectilignes dans le cristallin. 
Il faut pour cela supposer, pour chaque hauteur d’intersection avec 


(2) A. Brot, Ann. Soc. Scien. Brux., 74, 175-188, 1960. 
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la cornée par exemple, que l’indice reste constant mais est différent 
selon la hauteur. 

Ce mode de traitement ne suit évidemment pas la réalité car les 
rayons réels sont courbes dans l’hypothése d’une variation continue 
de l’indice. On voit de suite qu'il est toujours possible, si l’on opère 
de la sorte, de trouver un mode de variation de l’«indice » tel que 
toute loi d’aberration pourra être retrouvée. On voit aussi que 1’« in- 
dice » ainsi défini doit finalement diminuer quand la hauteur d’inter- 
section sur le cristallin augmente à partir de l’axe, de chaque côté. 
En effet, dans l’œil théorique, l’aberration est sensiblement sous-cor- 
rigée : les rayons sont trop réfractés. Pour qu'ils le soient moins, con- 
formément aux observations, il est indispensable que l’indice diminue 
et d’autant plus qu’on s’écarte de l’axe. 

b) Tenant compte de la courbure du rayon, on peut supposer, 
par exemple, que la loi de variation est linéaire pour l’indice ou le 
carré de l’indice (la différence est faible en pratique car le coefficient 
de proportionnalité est très petit). On prendra ici, comme en 18) dans 
le travail précédent 

= ne (1 + ay) 


no étant l’indice sur l’axe. Dans ce cas, le calcul numérique montre 
d’abord que pour obtenir une variation de l’aberration dans le sens 
d’une correction il faut prendre a=>0 c’est-à-dire que l’indice doit 
augmenter avec la distance à l’axe. On constate ensuite que la varia- 
tion de la hauteur d’intersection avec le plan focal est sensiblement 
constante quel que soit le rayon incident considéré. Or cela ne se peut 
car on obtiendrait ainsi deux «images» dans le plan focal. Il sera 
impossible de trouver une valeur de a qui rende compte des observa- 
tions. 
c) La loi de variation 
n= ne — by? 


qui donne dans le cristallin le rayon courbe 21) du travail précédent 
ne peut non plus être adoptée car on constate, par le calcul numérique, 
qu'elle conduit à une augmentation de l’aberration. 


3. Etude de l'application de la loi n?=n?+-by?. 

Il est utile de voir comment se font les calculs numériques. 

Le premier élément à connaître est la hauteur A, d’intersection à 
l’entrée du rayon choisi avec le cristallin. Cette valeur se déduit du 
calcul classique effectué pour ce rayon sur l’œil théorique normal. 
Mise à la place de y dans la relation qui donne l’indice on calcule la 
valeur me de ce dernier à l'entrée. On suppose que le cristallin possède 
cet indice et, par les formules classiques, on calcule la réfraction du 
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rayon incident. On possède de la sorte l’angle w, avec l’axe de la tan- 
gente à l’entrée au rayon courbe. On détermine alors la constante C 
de la réfraction C = n, sin à = ne cos We. 

On peut maintenant calculer l'équation du rayon. La forme géné- 
rale de cette équation est (I, (20)) 


(@: 
x + Cite = + nos Logio (y + 4/7? + (x — C?)/b) (1) 
Mb 

Nous connaissons M, C et b (que nous nous sommes donnés). 
Reste à déterminer la C‘*. Pour cela nous supposons que, pour x = 0, 
on ait y = yo. Il nous faut placer l’origine. Nous la placerons sur l’axe, 
a l’intersection de celui-ci avec la corde 2h, = 2yo. Par rapport au 
premier sommet du cristallin, l’origine se trouve au-delà, à une distance 

égale à la flèche qui correspond à cette corde, c’est-à-dire à 


fa = Rie — VR? — y? (2) 


a =0 


et la valeur de la constante s’écrit 


Cre = Logio (yo + / v2 -- (n; — C2?)/b) (3) 
My/b 
l’abscisse x1 de la deuxième surface du cristallin étant, en mm 
x1 = 4 — fh (4) 


Il est utile de transformer la quantité sous le radical dans (1) et 
dans (3). On a, en général, 
y? + (2 — C?) b = (n2 + by? — C2)/b = (n? — C2)/b = 
(n2 — n? cos? w)/b = n? sin? a/b (5) 


D’ou 


(Gis = ab 


Logio (Yo + Me sin Ge/b) (6) 


Mb 

C’est par cette relation que l’on calculera Ja valeur de la con- 
stante d’intégration. Dans le radical de (1) on calculera (nÿ—C?)/b. 
La valeur (5) nous servira encore plus loin. Les coefficients numériques 
de (1) sont maintenant connus. 

Il faut actuellement calculer les coordonnées de l’intersection du 
rayon courbe avec la deuxième surface du cristallin. Cela se fera par 
approximations successives. On se donnera par exemple une valeur 
approchée de l’ordonnée y; de sortie. Cette valeur pourra être celle 
que l’on a trouvée à l’aide du calcul effectué avec les éléments de l’œil 
théorique pour le même rayon incident. Avec cette valeur provisoire 
de ys on calcule d’abord la flèche f/2 correspondante sur le cristallin 
(Rec = 6mm) et l’on en déduit la valeur x11 de l’abscisse du point 


= sur le cercle 
Vea Vs DAMES oxi (7) 
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On introduit alors la valeur y = ys dans l'équation (1) et l’on 
calcule la valeur xj2 sur le rayon courbe. Si ys correspond à l’inter- 
section, on a x{1 = xio Dans le cas contraire on affectera ys d’une 
correction Ays que l’on pourra en général calculer avec avantage par 
la relation 


Ays = + (x1 — x12) tz We (8) 


Le signe à utiliser se détermine facilement à l’aide d’un croquis de 
la rencontre du cercle et du rayon. Le troisième essai donne en général 
des valeurs ys et xs suffisamment exactes. 

On calcule alors l’indice ns qui correspond à ys, puis par (5), la 
valeur ws. Celle-ci est l’angle que fait avec l’axe la tangente au rayon 
courbe à la sortie du cristallin. Pour achever de déterminer cette tan- 
gente, c’est-à-dire le rayon incident sur la seconde surface du cristallin, 
on calcule la valeur gs de l’abscisse par rapport au deuxième sommet 
du cristallin de l'intersection de cette tangente avec l’axe par la relation 


gs = + (vs/tgws + fle) (9) 
les signes étant à prendre d’après la géométrie du cas considéré. 


On est revenu au cas classique des rayons rectilignes. On trouve 
finalement A; dans le plan focal, la valeur que l’on cherche. 


4. Remarques. a) Si l’on détermine, en se servant de (6), la valeur 
de dy/dx, c’est-à-dire celle du coefficient angulaire de la tangente à la 
courbe au point (x, y), on trouve naturellement 


dyldx = + Vn®—C+y¥/C (10) 
Cette valeur est nulle pour 
r=4 (8 er 


Cela ne peut arriver que si n3} — C? est négatif. 

Le rayon ne possède un minimum dans le milieu cristallin que dans 
ce cas. Il est important de calculer alors, par (1), la valeur correspon- 
dante de x» sur le rayon. En pratique, deux cas peuvent se présenter : le 
minimum se trouve soit effectivement dans le cristallin soit en dehors. 
Dans ce dernier cas le coefficient angulaire de la tangente au rayon 
ne change pas de signe dans le cristallin; c’est le signe qui existe en 
Xo = 0. Ce cas se traite comme s’il n’y avait pas de minimum, et l’on 
applique directement la formule (1). Dans l’autre cas, (1) ne représente 
le rayon, dans le cristallin, qu'entre xo = 0 et x». Pour avoir l’autre 
partie de la courbe il faut introduire dans (1) le signe moins devant 
le radical. Cela provient de ce que, lors de l'intégration, on a posé 


VP + GB — C/b = 1— y 
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puis élevé au carré et enfin réintroduit cette seule valeur de 1—y après 
l'intégration alors qu’en fait t—y correspond à ++/y2+(n2—C2)/b. 

b. Dans le cas où l’on a (ng — C?)/b<0, on voit facilement que 
l’on a affaire à un minimum. En effet, puisque l’indice augmente quand 
y croît, la concavité du rayon est nécessairement tournée, en raison 
de la loi de la réfraction, du côté des y croissants. 

c. La valeur de ng — C2, we et ws peut devenir très petite. On n’ob- 
tiendra une précision suffisante qu’en prenant un nombre suffisant 
de décimales. Nous avons effectué nos calculs par logarithmes à 7 dé- 
cimales et c’est tout juste suffisant. Non pas tant dans le but d’obtenir 
des valeurs exactes que pour assurer la cohérence et l’exactitude des 
tendances de variation des résultats. 


5. Résultats des calculs. a) Les calculs ont été effectués pour l’in- 
dice nn et les 6 rayons méridiens utilisés dans le travail précédent. 
Les hauteurs d’intersection sur la cornée sont +1,75; +0.8; —0,2652 
(passant par le centre de l'iris); —0,6815 (perpendiculaire à la cornée); 
—2,0; —2,26. L’inclinaison commune de ces rayons sur l’axe optique 
est de 5°00’ 34’’,3. Pour les trois premiers nous avons pris Oy positif 
vers le haut; pour les trois autres vers le bas. Dans ces conditions, le 
calcul montre que le rayon commence par descendre vers l’axe dans le 
cristallin pour les n° 1, 4, 5 et 6 et monte pour les n° 2 et 3. 

b) Dans une première série de calculs nous avons pris b = by = 
{3,1319325} (2) = 0,00135498 qui donne pour le rayon 1 une valeur 
Ne = 1,42150 choisie arbitrairement. Les équations correspondantes 
des rayons dans le cristallin sont, dans l’ordre, pour les rayons définis 
ci-dessus 


x — 24,49393 = — {1,9489944} Logio(y + 4/2 — 3,133061) 
x + 29,99319 =  {1,9484441} Logio (y + 4/y? + 0,641078) 
x + 37,99608 =  {1,9474839} Logo (y + V7? + 7,204688) 
x — 49,36370 = — {1,9470177} Logio (y 1/y2 + 10,38023) 
x — 72,19765 = — {1,9451403} Logio(y + /y® + 23,10221) 
x — 74,56213 = — {1,9448165} Logio (y + a/y? + 25,28526) 


(1) Rappelons que {3, 1319325} représente le nombre dont le logarithme 
vaut 3, 1319325. 
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Pour le rayon | le minimum se trouve dans le cristallin. On a en 
effet xs = 3,57735 et Xm = 2.44341. Dans la première équation il faut 
mettre le signe + devant le radical avant le minimum et le signe — 
au-delà. 

c) Une seconde série de calculs a été faite en prenant pour b une 
valeur égale à la moitié de la précédente : b—b2 = }{4,8309025} = 
0,0006677449. Elle correspond, pour le rayon 1 à une valeur ne = 1,42072. 
Les équations deviennent 


x — 35,80502 = — {2,0992805} Logio (y + 4/72 — 3,123296) 
x + 53,94445 =  {2,0988943} Logio(y + y/y® + 2,171252) 
x + 72,59618 =  {2,0979989} Logio(y + \/y2 + 14,40908) 
x — 87,00042 = — {2,0975227} Logio (y + 4/7? + 20,89774) 
x — 116,95186 = — {2,0954818} Logio (y + 4/72 + 48,54393) 
x — 120,10006 = — {2,0951185} Logio (y + 4/92 + 53,43849) 


Ici aussi le premier rayon seul possède un minimum qui se trouve 
toutefois en-dehors du cristallin. On a en effet x; = 3,57832 et 
Xm = 4,72250. 

d) Remarquons que les valeurs de }; données dans le travail 
précédent correspondent à b = 0. 

e) Groupons en un tableau les valeurs trouvées pour /y. 


Rayon b, Différence b, | Différence | b— 

ot + : | 7 

| | 1,4606 0,0352 | 1,4248 0,0351 1,3897 

2 1,4937 0,0243 | 1,4694 0,0244 1,4450 

3 1,4680 0,0036 = 1,4644 0,0035 1,4609 

4 1,4519 | —0,0056 | 1,4575 | — 0,0054 1,4629 

5 14131 | —0,0317 | 14631 | —0,0318 1,4949 

6 | 1,4378 | —0,0346 1,4724 | —0,0345 1,5069 

| | 

EN ee eee 


6. Discussion des résultats. On voit d’abord (on reportera avec 
avantage les résultats sur un plan (Meornées At) que la courbe by se 
trouve a égale distance des courbes b; et b = 0 pour une méme valeur 
de Neornge : le déplacement des points de la courbe est proportionnel 
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a b. Cela permet de trouver facilement la courbe correspondante à une 
valeur utile quelconque de b. 

À partir d’une certaine valeur de b, il y a « correction » pour une 
ouverture déterminée de la cornée. La correction peut s'étendre au-delà 
de la zone optique de l’œil c’est-à-dire pour une pupille d’entrée plus 
grande que 4 mm en diamètre. Pour @ = 4 mm, la valeur optima de 
b paraît être un peu plus grande que by. 

Pour cette valeur il subsiste une sous-correction marquée pour 
le rayon | qui aborde l’œil le plus loin du côté temporal. II est intéres- 
sant de rapprocher ce fait du résultat de l’observation suivant lequel 
l’ouverture de l'iris n’est parfois pas centrée sur l’axe optique de l’œil 
mais est légèrement déplacée du côté nasal. 

Les trois courbes b = 0, bi, b2 se coupent sensiblement en un même 
point. Les coordonnées de ce point peuvent se déterminer par tâton- 
nements numériques. On peut en obtenir plus commodément la position 
par voie graphique en remplaçant les segments de courbe compris 
entre les points qui correspondent aux rayons 3 et 4, entre lesquels se 
trouve un point d’inflexion, par des segments de droite. On trouve 
ainsi pour l’intersection 

h 


cornée 


= 0,429 mm 


h, — 1,462 mm 


Ce point ne se trouve donc pas sur le rayon qui passe par le centre 
de l'iris. 

Il semble important de remarquer combien une faible différence 
d’indice entre l’axe et le bord du cristallin peut avoir d'effet sur les 
aberrations géométriques. On doit s’attendre à trouver, d’un individu 
à l’autre une forte dispersion dans les résultats des mesures de ces 
aberrations. Cette observation semble ne pas devoir être faite à l’occa- 
sion des mesures de l’aberration chromatique longitudinale. 


7. Conclusion. La loi de variation n? = ng + by? de l’indice dans 
le cristallin de l’œil théorique permet, malgré sa simplicité, d’aborder 
avec avantage l’étude de la représentation des aberrations géométriques 
observées sur l’œil réel en se servant de l’œil théorique. Il serait peut- 
être illusoire d’introduire un terme supplémentaire dans cette loi, au 
moins tant que les mesures d’aberrations effectuées sur l'œil réel 
n’auront pas abouti à des résultats concordants. De toute manière, 
il est nécessaire de tenir compte de la courbure des rayons dans le 
cristallin lors des calculs. 


oy! 
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a Bruxelles 


Premiére Section 


SCIENCES MATHEMATIQUES ET ASTRONOMIQUES 


Groupoides en relation associative 
et semigroupes mutuellement associatifs 


PAR 
Albert SADE 


1. Introduction. L'idée de considérer des groupoides en relation 
associative a été suggérée par l'étude des quasigroupes finis orthogo- 
naux, c.-a-d. satisfaisant 

VX), 2o0 APH BOR YS se KIS SS AVE I 
((8], p. 229). On sait que le produit de Mann E(+) ({“], p. 79, N° 8) de 
deux quasigroupes orthogonaux E (.) et E(x), 
xey = (x.y) X y 
est encore orthogonal à chacun d'eux. Dès lors, si un produit comme 
(x.y) X z pouvait être remplacé par x.(y X 2), il en résulterait une grande 
souplesse dans les calculs où intervient le produit à droite. On est ainsi 
conduit à dire que deux groupoides G = E(.) et H = E(x), sur le 
méme ensemble E, sont en relation associative si 
wx,y,zEE, (x.y)xz=x.(yxz), 
ce qui sera noté G (AM?) H. 

Cette relation demosienne s’exprime également par la permutabilité 

des translations 


IS As = NT TE 
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Si l’on a en plus 
@ ¥y).2=x xX &.2Z), c-a-d.. AZ) CG: 
alors on dira que G et H satisfont à l’associativité mutuelle, notée G 
(S HN) H. 
On trouve de telles relations, par exemple, dans le cas des rin- 
goides ([°1, p. 47). Les définitions et notations sont celles de [*] et de [5]. 
ll ne sera question que de groupoides homogènes. 


2. Exemples. (1) Sur un anneau quelconque, A, les groupoides 
X.Y =f(X) y, 
eX Ce a), 
ou fet g sont deux fonctions arbitraires définies sur A, à valeurs dans A, 
sont (%.%). 
(ii) Le semigroupe de Thierrin, xy = y ((4], p. 153, N° 2, équ. 9) 
est (#.A) avec tout groupoide, défini sur le même ensemble. 


(iii) Soit G = E(.) un semigroupe abélien d’ordre n sur un ensemble 
E de n éléments. Pour que les groupoides 


P= EBs), eX y=] a. x y 
O) = IE 352 Issn! 


a, b, x, y EE. avec m, p, q, r arbitraires, non nuls, fixes € Z/n, soient 
(AA) il faut et il suffit que a.b? = a".b dans G et que m = mr, q = qp, 
dans Z/n. 


3. On vérifie sans peine que 

(i) Si deux groupoides quelconques P = E(.) et DEC) sur un 
même ensemble E sont (#-V) et s’il existe un élément u € E neutre 
à droite dans l’un et à gauche dans l’autre, alors P et Q coincident. En 
particulier, deux loops (#7) avec la même unité, sur le même ensemble, 
coincident. 

(ii) SiJ C E satisfait y x € E à la condition J.x = J (resp. x.J = J) 
sur P, alors J satisfait à la condition J x x =J (resp. x X J = J) sur Q; 
P et Q étant (2H). 

(iii) Pour qu’un groupoide G = E(.) et son conjoint C = E(x) 
((4], p. 155, ii), définis sur un même ensemble quelconque E. satis- 
fassent G (Z.7) C, il faut et il suffit que G obéisse à la loi de Hosszu 
ayaa ty. (xz) (7) p: 205). 

(iv) Pour qu’un quasigroupe a gauche G = E(.) et son réciproque 
Q = E(x) soient (7.7) il faut et il suffit que G soit un semigroupe. 


4. Construction. Si, sur un ensemble quelconque E, P = BC rest 


sh 


un semigroupe avec unité a droite en relation associative avec un 
groupoide Q = E(x) sur le même ensemble, P (47) Q, alors . 
(i) la loi de Q est x X y= x.f(y) où f est une fonction arbitraire, 
définie sur E et prenant ses valeurs dans E, et réciproquement. 
(ii) si, de plus, Q est cancellable a gauche, alors Q est l’image de P 
par Visotopie T = {2 = 1, n = y + g(y), € = 1}, dont les composantes 
£ et C sont réduites à l'identique et où g est la fonction inverse de /. 


ty 
— 


Preuve. Par hypothèse, (x.y) X z = x.ly X 
(i) Soit u l'unité à droite de P; alors 


CDR EE OU Ce 


or u X z ne dépend que de z puisque u est invariable; en posant 


HSE), 
on à 
RS ee ACIER 
où fest une fonction arbitraire de z. Réciproquement, si x X z = x.f(2), 


puisque P est un semigroupe, on aura 


(xy) % 2S Cop he = Er JE. 
D’autre part 

(vie 2) = (ye) ee 
(i) Si Q = E(X) est cancellable à gauche, 


f(z) =e) = 2S kh 2s Ses 


et f est une permutation de E. Elle a done une fonction inverse g(z) 
et Q est le produit de P par l’isotopie de composantes = 1, y = 
(y > g(y)), = 1. On a un énoncé analogue si Q est un semigroupe 
avec unité a gauche. 


5. (i) Deux groupoides P = E(.) et Q = E(X), construits sur le 
même ensemble E, sont en relation associative, P (#.o7) Q, si T = {E=1, 
n—=C—x->x X m} où m est un élément fixé quelconque € E, est 
une autotopie de P, quel que soit m. 

(ii) Un groupoide P = E(.) admet un groupoide Q = E(x) en 
relation associative. P (4.47) Q, si le groupe d’autotopie de P a un sous- 
groupe de la forme (— = 1, » <= (); 

(in) Pour qu'un quasigroupe P = E(.), construit sur un ensemble 
quelconque E, admette un groupoide Q = E(x) en relation associative, 
P(A.S) Q, il faut et il suffit que P soit isotope d’un groupe quelconque G. 
Alors Q est un quasigroupe et P et Q sont respectivement isotopes de G, 
avec les composantes (&, n, n) et (n, À, 7). 
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Preuve (i) Si T est une autotopie de P et si x.y = z, on aura 
Om) 2 Gs 
c’est-à-dire 
OM BOS Ww = Cenc iid 

Cette égalité étant vérifiée quels que soient x, y. m, le groupoide 
Q est en relation associative avec P. 

(ii) Les autotopies de la forme £ = 1, n = €. forment évidemment 
un groupe, celui-ci étant isomorphe au groupe constitué par les com- 
posantes 7. 

(iii) La proposition est une conséquence du théorème ([$], p. 632 
N° 7.2; [6] p. 304). On peut toutefois prouver directement que, s’il 
existe un élément u, commun à deux groupoides (#V) P = E (.) et 
Q = E(x), non singulier à droite [>], dans l’un et à gauche dans l’autre, 
alors P et Q sont isotopes d’un même semigroupe avec élément neutre. 


6. Produit cartésien. Soit (Ei), i € 1, une famille d’ensembles, P: 
et Q; deux groupoides définis sur E; et wi EI, Pi (AV) Qi. Alors, si 
P est le produit cartésien des P; et Q celui des Qi, définis sur le produit 
cartésien E des Ei, on a P (Z.:27) Q. 


7. Produit selon G(o). Dans [3], N° 4, p. 230, on a défini le produit 
S = E(#) de deux groupoides quelconques P = E(.) et Q = E(x), 
selon un groupoide G = E(o) par la relation 
xey = (x.y) O(xxy) , Vx, VEE. 
Si P et Q sont (Z.7) sur E, et si S = E(*) est leur produit selon 
G = E(0); (i) quand P (4.7) Q et quand P est un semigroupe distributif 
à gauche par rapport a G, alors P (#.A) S; (11). SiQ (ZA) P et si Qest 
un semigroupe distributif 4 gauche par rapport a G, alors Q (42%) 8; 
(iii) Si Q (AA) P et si P est un semigroupe distributif a droite par rap- 
port à G, alors S (ZA) P; (iv) Si P (ZH) Q et si Q est un semigroupe 
distributif à droite par rapport à G, alors S (4.7) Q. Démontrons seule- 
mer t (i). 
Par définition 
(x.y)#z = [(x.y).z] OL@.y) xz] 


= [x.(y.z)] O[x.(7 xz)] 


x.(y#z) = x.[(y.z) O(v Xz). 
A cause de la distributivité de P les derniers membres sont égaux. 
L’égalité des premiers membres exprime que P (2D) S. 


8. Produit à droite. Appelons produit à droite (Généralisation du 
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produit de Mann) de deux groupoides P = E(.) et Q = E(x) le grou- 
poide T = E (#) défini par 
xæy = (x.y) X y 3x, yEE. 
Un groupoide est diagonal si l'équation x.x =a a une solution 
unique en x, Ya. | 
Si P (AA) Q, (i) pour que P (Z.2) T il faut et il suffit, si Q est 
diagonal, que P soit associatif; (ii) pour que Q (ZA) T il faut et il 
suffit que Q (AA) P; (iii) Si Q est self-distributif, alors T (AA) Q. 
Par définition on a 
Trt] fe eh) Fg Eee a Bcd hoe 
et comme P (#.A) Q 
AT (EE). 
Ensuite 
B =x. Ora) = RIM IE 
Pour tout 1 EE, dz, z X z—1, puisque Q est diagonal; donc 
P(Z.97) T= A =B= TT XP. 


et P est un semigroupe. La preuve est analogue pour les autres parties. 


9, Homomorphisme. Soit T une application de l’ensemble E sur 
l’ensemble E’ = Er telle que l’image par T du groupoide P = E(.) 
soit P* =, PR et celle de Ge Hic), O° = OR Alors 


P (PA) Q = P' (PA)Q’. 


10. Zsotopie. Soit P’ l’image de P = E(.) par l’isotopie de compo- 
santes (&, C, ©) et Q’ celle de Q = E(x) par l’isotopie (&, n, 0). Alors 
P (RA) Q = P’ (#A) Q’ 
om 


où &, n, © sont des permutations arbitraires de E. 


11. Transitivité. La propriété d’être en relation associative est 
transitive car, soient P = E(.), Q = E(x), S = E(*) trois groupoides 
(homogènes) avec P (AA) Q et Q (AA) S; donc 


(x.y)Xz = x.(yXz) ; (xxy)az = x X(yxz). 
Par suite 
[x.(y X z)]#t = [(x.y) x z]at = (x.y) X (zxt) 
= x.[y X (zxt)] = x.[(y x z)*Z]. 
Or V_ a, By, 2, y X z=u, donc 


(XS. ust — x. (et): 
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12. Associativité mutuelle. La plupart des propositions précé- 
dentes restent valables pour l’(% -@). Si deux groupoides (4.2) sont 
abéliens il est évident qu'ils sont mutuellement associatifs. Pour que 
deux groupoides soient en (.7 -@) il faut et il suffit qu'ils soient asso- 
ciatifs. 

On trouve encore facilement que si un groupoide cancellable d’un 
côté et son conjoint sont en (.7% .#) ils coincident. 

Pour qu’un quasigroupe à gauche G = E(.) et son réciproque 
Q = E(x) soient (2 .#), il faut et il suffit que G soit un semigroupe 
satisfaisant à l’équation de transitivité a.b = (a.c).(b.c), ([{], p. 203). 
G posséde alors une unité a droite unique et tous ses éléments sont 
unipotents. Si, de plus, G est cancellable, à gauche, alors il est en outre 
commutatif. 

Si P= E(.) et Q= E(X) sont deux semigroupes en (.7 .#1, 
G = E(o) un groupoide entropique quelconque, sur le même ensemble, 
et si P et Q sont distributifs des deux côtés par rapport a G, alors 
(i) leur produit H = E(*) selon G est encore distributif bilatère par 
rapport à G, (ii) H est en (WV -@) avec P et Q. 

La relation d’associativité mutuelle est une équivalence par rapport 
à l’ensemble S de tous les semigroupes construits sur un ensemble 
commun E. La partition modulo (7 .#) de S est invariante par toute 
permutation de E. 
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Un modéle mathématique de systeme a dissipation 
du type «structurel » 
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F. BUCKENS 
SOMMAIRE 


Le modèle mathématique d’un système oscillant à un degré de liberté, 
à pertes « structurelles », est généralement représenté par une équation linéaire 
«approchée » du second ordre dont le coefficient d'amortissement est choisi 
de façon à donner, en régime périodique, une dissipation cyclique propor- 
tionnelle au carré de l’amplitude et indépendante de la fréquence imposée. 
On sait que le sens physique de pareil modèle linéaire où le coefficient d’amor- 
tissement devrait être inversement proportionnel à la fréquence imposée, 
s’explique par la considération d’un système à amortissement « héréditaire », 
correspondant en réalité à une équation d'ordre plus élevé, mais dont la 
linéarité permet, grâce à une transformation de Laplace, de réduire la dis- 
cussion à celle d’une fonction de transfert ou, si l’on veut, à celle d’un module 
complexe. Si le schéma est toujours possible à l’égard d’oscillations harmoni- 
ques imposées en prescrivant que la dissipation cyclique est une fonction 
donnée de la fréquence, il perd son sens physique pour les oscillations tran- 
sitoires lorsqu’on suppose que cette dissipation est complètement indépen- 
dante de la fréquence. 

La recherche directe d’une force d’amortissement satisfaisant à la con- 
dition des pertes cycliques à l'égard des fréquences sous la forme d’une fonction 
généralisée au sens de Temple-Lighthill, conduit ici, sans introduire de termes 
de mémoire, à une expression simple, mais effectivement non linéaire de 
l’équation différentielle. La reconstitution des cycles d’hystérésis classiques 
par l’introduction d’un type assez particulier de mémoire « courte » conduit 
à une forme plus générale dont la connaissance est précieuse pour l’analyse 
ultérieure de modèles mathématiques corrects sur machine analogique. 

En guise d’application, l’étude des oscillations naturelles amorties, dans 
le cas particulier d’une forme exacte mais à linéarité « fragmentaire» et 
symétrique du modèle mathématique, sous la forme : mx + k(x + Bxlx/x) 
= 0, donne des résultats simples mais exacts qui font apparaître, par rapport 
à ceux fournis par l'emploi de l'équation linéaire « approchée », de grosses 
différences, tout au moins en ce qui concerne la variation de la pseudo-pulsa- 
tion en fonction du coefficient d’amortissement 8. 


SUMMARY 
The mathematical model of an oscillating system with one degree of 


freedom, whose energy dissipation under forced oscillations is of the structu- 
ral type, is usually given by a second order linear equation, the damping 
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coefficient is then defined by the condition that cyclic dissipation should be 
proportional to the square of the amplitude and independent of the frequency 
(internal or structural damping). The physical meaning of such a linear model, 
where the damping coefficient should be inversely proportional to the imposed 
frequency, receives its natural interpretation by introducing an « hereditary » 
damping, and thus higher derivatives in the equation, and by considering 
the transfer function after a Laplace transformation. 

Such a scheme allows simulating any cyclic dissipation as far as harmonic 
oscillations are concerned, but runs into difficulties for transient motions, 
when the dissipation is assumed of the «structural type», i.e. completely 
independent of frequency. 

The search for a damping « force», with cyclic dissipation independent 
of frequency under the form of a generalised function (in the sense of Temple- 
Lighthill), leads to a simple but non-linear expression of the differential equation. 
The further introduction of «memory» terms of a special type allows re- 
shaping the solution in accordance to the classical hysteretical cycles, and 
this lead to a more general form of the model. Its knowledge is valuable, for 
it makes possible analyzing correct models on analog computers. 

As an application, the damping of natural oscillations is studied in the 
particular but exact case of a «fragmentary» and symmetrical linearity of the 
mathematical model in its simplest form: mx + k(x + Bx\x/x}) = 0. The 
solution is straightforward and shows large differences with results obtained 
from the «approximate» linear model in the behaviour of the pseudo-fre- 
quency as a function of the damping coefficient 8. 


1. DISSIPATION STRUCTURALE 


La littérature technique consacrée à la dissipation de l’énergie 
vibratoire par hystérésis fait généralement intervenir un module com- 
plexe d’élasticité dont le rôle restrictif ne semble pas toujours avoir 
été bien compris; une utilisation inconsidérée de ce module peut con- 
duire à des erreurs importantes comme on peut le voir par des déve- 
loppements mathématiques assez élémentaires. 

La linéarité complète de l’équation du second ordre du mouve- 
ment d’un oscillateur soumis à une force extérieure P cos (wt), variant 
sinusoidalement avec le temps, ne s'obtient que sous la forme : 


mx + cx + kx = P cos (wi) (1) 


où m, c et k représentent l’inertie, le coefficient d’amortissement et le 


module de raideur. 
La solution de régime : 


x = a cos (wt — 9) (2) 


comporte une amplitude relative a/ao, (où ao = P/k), dont l'expression 
ESLÉE 


a/ao — al + 02)? ae 46292] 1/2 (3) 


a 


et un retard de phase donné par : 
tan = 2Cp/(1 — p?) (4) 


Dans ces expressions p = w/wo (où wo = Vk/m), et © = c/2mao. 
L’énergie E dissipée par cycle est, avec T = 2x/o : 


T T 
E = Ap P cos wt.xdt = [ cx? dt (5) 
0 0 


E = xPa sin (®) = rca. (6) 


d’où : 


Cette dissipation est ici proportionnelle à la fréquence et au carré 
de l’amplitude. 

Or, des expériences déjà anciennes conduites par Kimball et Lovell[i] 
ont montré que pour de nombreux matériaux, cette dissipation cyclique 
est bien proportionnelle au carré de l’amplitude, mais pratiquement 
indépendante de la fréquence, et ceci pour une très large bande de fré- 
quences. On peut donc parler d’une dissipation par hystérésis. C'est en 
effet ce que l’on retrouverait en admettant que l’aire entourée par la 
courbe tracée dans le diagramme tension déformation est indépendante 
du temps mis à parcourir le cycle. 

Des expériences plus modernes (Lazan [?]) confirment ces résultats 
en restreignant quelque peu le domaine des fréquences et des ampli- 
tudes. 


2. LE MODÈLE CLASSIQUE, LINÉAIRE MAIS APPROXIMATIF 


En invoquant la petitesse de l'effet, l'habitude s’est installée de 
parler de coefficient d'amortissement inversement proportionnel à la 
fréquence imposée : 


CPs (6) 


de manière à conserver une solution pratiquement sinusoïdale qui 
donne à la fréquence considérée, une dissipation cyclique indépendante 
de tr: 


E = tha? (7) 


Comme la constante / a les mémes dimensions que k, il est d’usage 


d'introduire un paramètre 1 sans dimension, en posant : h = yk. On 
a alors 


6 = h/2moow = nwo/20 = y/20 
ce qui donne le modéle « approché » : 


mxX + (h/w)X + kx = P cos (wf) (8) 
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ét au heu de (3) et de:(4) « 
alao = [(1 — 92)? + n?} 1/2 (9) 


tan @ = »/(1 — p?) (10) 


Il est clair que la forme (8) de l’équation n’est pas satisfaisante 
et qu'il convient de reconsidérer le problème sous l’angle de la trans- 
formée de Fourier si l’on conserve l'hypothèse de la linéarité. 

En admettant alors que le degré de l’équation du mouvement 
puisse être supérieur à deux (système à plusieurs degrés de liberté), ce 
sont les dérivées d’ordres impairs qui contribueront à former la partie 
imaginaire Ay du module complexe k* = k (i + jy). Celle-ci doit être 
une fonction impaire de w, prenant donc ici dans le cas actuel, la forme 
kno/|o!. 


On peut alors montrer (Fraeys de Veubeke [?]) que l’on aboutit 
à une absurdité physique en régime transitoire du fait que l'effort 
interne dépendrait aussi bien de l’histoire future que de l’histoire passée 
du mouvement. ” 

On peut échapper à cette impasse en supposant que la condition 
sur l'indépendance de l’amortissement cyclique vis-à-vis de la fréquence 
ne vaut que dans une bande plus ou moins large; Biot [°] donne une 
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solution approchée valable dans un domaine ne contenant pas le voi- 
sinage de wo = 0. 

Nous cherchons ici une solution qui maintienne intégralement le 
caractère d’indépendance de la dissipation cyclique vis-à-vis de toutes 
les fréquences, en laissant donc tomber nécessairement la condition de 
linéarité de l’équation. 

Nous admettrons que cette dissipation ne dépend que de l’am- 
plitude a de l’oscillation, que nous supposerons d'abord complètement 
symétrique. Il n’est pas nécessaire de spécifier dès le début la forme de 
cette fonction F(a). Remarquons d'ailleurs que certains résultats expéri- 
mentaux indiqueraient, pour les poutres composites, une dissipation 
proportionnelle au cube de l’amplitude de déformation [°], et d’autre 
part, on sait que dans des phénomènes d’hystérésis magnétiques, (fig. 1), 
Steinmetz adopte une perte par cycle indépendante de w et proportion- 
nelle à al8. 


3. DÉTERMINATION D'UN MODÈLE EXACT, A DISSIPATION STRUCTURALE, 
SANS MÉMOIRE 


Reprenons donc le problème de plus haut et, en admettant que la 
dissipation cyclique soit indépendante de la fréquence et donc seule- 
ment fonction de l’amplitude des oscillations entretenves par une force 
extérieure P cos wf, demandons-nous quelle doit être l'expression de 
la force interne suscitée par le mouvement, sans chercher à séparer dès 
l’abord la partie conservative de la partie dissipative. En désignant cette 
force par — & (x, X), nous avons |’équation : 


mx + @ (x, *) = P cos wrt (11) 


et la condition : 
= Hig 
E = Î P cos wt.Xdt = f A(x, x). xat = Fa) (12) 
0 0 


avec F (a) — 0, pour toute amplitude a et toute pulsation imposée w. 

En effet, au lieu de P et w, nous prendrons a et «© comme para- 
mètres indépendants, a étant une fonction de P et de w, ainsi que des 
coefficients de raideur, amortissement et inertie du système. Ces der- 
niers coefficients sont pris fixes et dès lors, nous pouvons écrire : 


x(t) = a. ET, p, a) (13) 


où + = t/T, p = w/wo et « = a/ao sont des paramètres sans dimensions. 
Le choix des temps et longueurs caractéristiques 27/wo et ao (éventuels) 
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pourra se faire en prenant par exemple : 
wo? = © 3 (0,0)/m, ou [a (0,0)/S° (0,0)}>, ou @ “(0,0)/ @ ” (0,0), etc. 
(14) 


" 


a, =%S,,0,0)/2/(0,0) ou 2”, (0,0)/2; (0,0), etc. 
we “ TT Go 
le choix entre ces différentes expressions étant indifférent pourvu qu'il 
soit significatif (différent de 0 et de co). 
Dès lors, les expressions (11) et (12) peuvent s’écrire (en posant 
Déjor = &) : 


(p/2m)2ak + 2 (ak, ab w/2m)/k + aË = (PJk) cos (2x7) (15) 
af . L = = 
Bia =P f cos (2m T)Edt = Hi @ (aë,aËw/27) Edt = F(a)/a (16) 
0 0 


L’équation (15) montre que d’une manière générale & dépend 
bien de ©. D'autre part, d’un point de vue purement dimensionnel, 
nous devons avoir, en dérivant par rapport à @ : 


L : : 
| S: (a, a&w/2m) &2 dr = 0 (17) 
0 


En fait, l'intégrale comprenant #7 ne devrait être étendue à l’inter- 
valle (0,1) qu’en exceptant les points où £ — 0, pour lesquels cette dérivée 
pourrait ne pas exister; la présence du facteur £2 permet en fait d’intro- 
duire des fonctions généralisées de la forme S(Ë), 3'(Ë), ... (fonction 
de Dirac et dérivées) dans l'expression de @z, sans changer l’intégrale 
complète (17). 

Une nouvelle dérivation de (17) par rapport à w donne : 


1 5 
[ oi É3dr =0 (18) 
0 


etc., de sorte que l’on a nécessairement 2% Dino =... = 0 pour tout 
x et pour £ + 0, si 2(x,x), et donc aussi 2% (x,*) appartient à la classe 
des fonctions généralisées au sens de Temple-Lighthill, [6], [“], (voir 
Annexe 1). 

On est ainsi conduit a écrire l’expression : 


Bi (x, à) = fi (x, X) 8%) + fax, À) SC) + … (19) 


La détermination des fonctions fi (x,x), f2(x,*), ... ne peut résulter, 
dans notre cas, que de la condition (12), jointe aux conditions de symétrie 


postulées ici. | | 
L'intégration de (19) opérée dans les deux sens de x à partir de 
x = 0, donne, en introduisant g(x) comme «constante » d'intégration 


63 


par rapport a x, et en posant : 


4h(x, X) =[f, (À, À) — fi X)] (20) 


G(x) = g(x) + A(x,0+) — h(x,0 —) (21) 
l’expression (v. Annexe 2) : 
@(x, X) = [A(x, 0) + h(— x, 0)].X/|X| + [A(x, 0) — h(— x, 0)]x/|x| 


+ [fo(x,0-+) u(x) — fe (x,0—) u (— X)] 8(X) + GO) (22) 


Il reste a vérifier la condition fondamentale (12). Par substitution 
llerestes: 


gite a 
f [A(x,0) + h (— x,0)] |x| dt = 4 if [h(x,0) + h (— x,0)] dx = F(a) (23) 
0 0 


à satisfaire pour tout a> 0. On a donc : 

h(a0) + h(— a,0) = F'(a)/4 (24) 
pour tout a > 0, et donc : 
4 [h(x,0) + h(— x,0)] = 


RGO) + LC x,0) — fe, 0) — £5, (= 2,0) = Fx) (25) 


Comme nous avons supposé, dans ce qui précède, que le système 
est entièrement symétrique en ce qui concerne ses oscillations positives 
et négatives, nous ne pouvons garder ici que les termes qui changent de 
signe lorsqu'on change simultanément les signes de x et de x. 

Ceci élimine le second terme de (22), qui contient une fonction 
symétrique en x; le troisième terme ne peut être conservé que si 
Jax,0+) = fo(—x, 0—) = f(x), pour tout x. Enfin, la fonction G(x) 
devra elle-même être impaire g(x) = — g(— x). 

Le modèle mathématique cherché pour le système symétrique à 
dissipation structurale, sans termes de mémoire, est donc représenté 
par l'équation : 


mX + F’(\x|) X/4 |X| + f(x) u(X) — f(x) u (— XJR) + g(x) =P coswt (26) 
dans laquelle la force de rappel g(x) est impaire. 


4. TERMES DISSIPATIFS ET CONSERVATIFS. EXTENSION AUX SYSTÈMES 
DISS YMETRIQUES 


On voit qu’au second membre de (22) seul le premier terme (25) 
donne une dissipation cyclique; on peut le considérer comme le terme 
caractéristique de dissipation structurale. 


Le second terme fait apparaitre dans le crochet une fonction im- 
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paire en x; ce terme correspondrait à une résistance tantôt positive, 
tantôt négative à dissipation totale nulle: en fait, sur chacun des inter- 
Ales Det 0:ona: 


+a 0 
( [ + [ ) to) n= x0) as] = He@ = Hea (21) 
0 +a 


Du point de vue physique, il suffirait que : 
h(x,0) — h(— x,0)< F’(\x\)/4 (28) 


pour que la résistance provenant des deux premiers termes ensemble 
ne devienne jamais négative. Or, d’aprés (20) et (25), ou encore d’aprés 
(k), Annexe 1, cette condition se ramène a: 


h(x,0) > 0 (29) 
pour tout x, positif ou négatif, soit, d’aprés (28) : 
f(,0) — f,,,(%.0) > 0 (30) 


Bice 


La figure (2) donne les courbes de F’(|x|) et de Hi) = 1620) 
— h(—x, 0) (qui a été supposée continue), ainsi que celles de leur somme. 

Notons que si le second terme de (22) ne donne aucune dissipation 
cyclique, il ne peut être utilisé dans le modèle symétrique. 

Le troisième terme représente une force d’un type assez Curieux, 
nulle tant que l’oscillateur est en mouvement (x * 0), prenant une valeur 
instantanée, éventuellement très grande lorsque la vitesse s’annule 
(x = 0), et de plus, changeant de valeur d’après que la vitesse rejoint 
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ou quitte la valeur 0 par valeurs positives ou négatives, pour autant que 
f(x, 0+). fe(x, 0 —) > 0. Il s’agit donc d'une force de « frottement 
statique », dont la valeur pourrait dépendre non seulement de l'écart 
x = +a, mais encore du sens du mouvement x = 0 +. Le coefficient 
de frottement statique ne serait symétrique que si f2 (x, 0 +) = f2(x,0 —), 
soit pour le système symétrique (26): f(x) = /(— x), fonction paire. 

Cette résistance au mouvement disparaît dès que x 0, ce qui 
explique qu’elle ne dissipe aucune énergie, tout en étant évidemment 
capable de modifier le mouvement. 

Enfin, le dernier terme n'est pas dissipatif si G(x) est une fonction 
univoque de x. 

On peut songer à étendre les résultats précédents aux systèmes 
dissymétriques, en mettant par exemple la condition (12) sous la forme : 

E = F(a’) + F(a”) + 4 [S(a’) + S(-a’) — S(a”) — S(—a") (31) 
où a’ et — a’’ désignent les amplitudes positive et négative. 

En définitive, si nous laissons de côté des forces de nature assez 
particulière n’intervenant qu'à l’arrét (* = 0), telle celle représentée 
par le troisième terme de (22), le modèle mathématique d’un système 
oscillant sans mémoire, à dissipation cyclique caractérisée par (31), 
répondrait à l'équation 

; EL nt Len Neer , 
mx + Pus )+ IT) x/2|x| + G(x) = P cos wt (32) 
dans laquelle G peut être dissymétrique. 

Si l’on admet que la seule force responsable de la dissipation 
cyclique est «fragmentairement » linéaire, c’est-à-dire qu'elle est 
linéaire dans les intervalles (x > 0, x > 0), (x > 0, x < 0), (x < 0, 
X <0) et (x < 0, * > 0), une solution pourra s’obtenir par raccorde- 
ment de solutions élémentaires construites dans chacun de ces inter- 
valles si les troisième et quatrième termes présentent le même caractère 
de linéarité fragmentaire, c’est-à-dire que dans l’équation : 
ona: 


MX + [K’u(x) + K"u (— x)] X |x/X| + H(x) x /|x| + G(x) = P cos wt 
H(x) = Hox/|x| + Hix 


G(x) = Gox/|x| + Gix 
etvalors : 


E = K’a” + Ka" + Hi(a’ — a”) (33) 


Les solutions partielles sont sinusoidales avec chacune deux con- 
stantes arbitraires a fixer par les conditions de raccordement (en x et x) 
entre les quatre intervalles précédents (i = 1, 2, 3, 4) 


xi = coi + [P/m(w — wi)] cos wt +Aicos(wit) + Bisin (cit) (34) 
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avec co. = — (Go + Ho)/(G1 + H1 + K’, 


coz = — (Go — Ho)/(Gi — Hi — K’) , 


cos = — (Go — Ho)/(Gi + Hi — K”) , 


I 


— (Go + Ho)/(Gi — Hi + K”), 


Co4 


5. OSCILLATIONS LIBRES D'UN OSCILLATEUR 


L'étude de l’amortissement de l’osci 


Os = G; + Hi + K’ 


mon = G1 — Hi — K’ 
(35) 


9 
mas = Gi AK" 


2 
MmM@O4 = Gi — Hy + Ke 


SYMÉTRIQUE 


llation naturelle à partir de 


conditions données (xo, Xo) se fait très simplement à partir du schéma 


précédent, en posant P = 0. 


Limitons-nous ici au cas du modèle sous sa forme la plus dépouillée, 
po (55 — Hp — Hy = 0 Gr het Ke REF 


PES 
SS 
de CS 
D 
Si 
| 
’ | 
T/4 4 
Fig. 3 
L'énergie dissipée par cycle serait ici; d’après (33) : 
E = F(a) = 28ka? 
L’équation est donc alors : 
mx + k(x + Bx|x/x|) = 0 (36) 


Nous avons donc coz = 0, et : 


oF = oF =o2%=k(1 + 8)/m= or +68), [’ =2nr/o’ (37) 


2 


Wis wo? Bis"? 2 k(1 B)/m = or = 6) à F7 = 2Tr/o” (38) 
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Partons de to = 0, avec l’amplitude Ai = a); ona successivement, 
pour les intervalles : 
Gee 7S hy = TS 
{1 = ar! <—f—t—(1 = 1/4 
to =(T’ + T)/4 <t<t3 = (2T" + T’/4 
ta = COTE TA EI Et = (T+ 77/2 
les solutions respectives : 


x(t) = ao cos w”1 

x(t) = — (a o)”/@’) sin w’ (t — f1) 

x(t) = — (ao w"/@’) cos ©” (t — te) (39) 
x(t) = (ao w"?/’2) sin w(t — fz) 


pour lesquelles : 


x(t) = — ao o” 
x(t2) = — aow"/w’ = — ar 
X(t3) = ao w"?/@’ = aaa” 
x(t4) = ao w"2/@’? = a2 0"/@’ = aa 

La figure 3 représente cette oscillation amortie, composée de 
tronçons de quarts d’ondes sinusoidales, le coefficient de décroissance 
des amplitudes au bout d’une pseudo-période T = (T’ + T’’)/2 vaut 
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donc w’/w’? = (1 — B)/(1 + 8), et la pseudo-période étant : 
T = (x/o0) [( + 8) + (1 — B)-1/2] (40) 
la pseudo-pulsation d'amortissement vaut : 
@a = 2n/T = 20/[(1 + B)-1/2 + (1 — B)-1/7] (41) 
Si le coefficient d'amortissement B est petit (en tout cas < 1), on peut 
employer le développement en série : 
Ga = wo [1 — 382/8 — 1784/128...] (42) 


La fig. 4 représente la méme oscillation dans le plan de phase 
(x, X) où la «trajectoire » s’obtient en raccordant des quarts d’ellipses. 

Il est intéressant de comparer le résultat (41) avec celui que Collar [8] 
obtient en admettant la validité de l’équation linéaire (1), avec P = 0 
et en y substituant d’après (6) € = h/wa = nk/oa, Où wa est précisément 
la pseudo-pulsation à calculer. Il trouve ainsi : 


(oo) cn, / (i + 4/1 — »)/2 (43) 
dont le développement est : 
(wa) = wo[l — 72/8 — 5y4/128...] (44) 
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La comparaison doit se baser sur l’égalité des dissipations cycli- 
ques (7) et (12), en régime périodique; on en déduit mh = 28k d’où 


= rn/2 (45) 
ce qui fournit l’expression exacte (43) : 
Ga = @o [1 — 3(r/2)? n°?/8 — 17(x/2)4 4/128...] (46) 


Fig. 4 


et ceci donne pour la pseudo-pulsation «4, une décroissance avec le 
coefficient n beaucoup plus rapide que dans la formule de Collar. 

D'ailleurs, la formule (42) donne oq = 0 pour B = rn/2 = 1 
c’est-à-dire pour une valeur de 7 qui ne rend pas encore Vexpression 
(44) imaginaire. Ceci s’explique si l’on note que d’après (38), on a 
dans ce cas w’’ = 0 de sorte que x(t) = ao = constante dans l’inter- 
valle infini 0 < t << t =—T’’/4 = w: le système conserve écart x = ao 
car la résistance au mouvement équilibre exactement la force de rappel, 
kao. 

D'autre part, on note que si ® < 1, mais proche de l’unité, le 
coefficient de décroissance (1 — 8)/(1 + 8) provoque une chute pres- 
que totale au bout d’une seconde période, d’ailleurs longue. Si B « 1, 
on vérifie que cette loi de décroissance donne, jusqu’au second ordre, 
le même résultat que celle fournie par la méthode approchée de Collar; 
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celle-ci donne un facteur de décroissance 


e-eL/2m = e-mn/V(141—n)/2 = | — yr Jo + 4/1 — y)/2 


+ nn +4/1 — x) — … 
alors que le coefficient exact vaut : 
(1 —B)/ +8) = 1 — 28 + 26? +... (47) 


La comparaison donne : 


B = (y7/2)/ J + 4/1 — 4)/2 = (yr/2)(1 + 42/8 +...) (48) 


et donne une coïncidence jusqu’au second ordre avec l’expression (46), 
obtenue par l’équivalence des énergies dissipées dans des conditions 
différentes, en oscillation forcée périodique. 

Notons cependant qu’en bonne logique,la méthode de Collar devrait 
faire intervenir dans le coefficient d’amortissement c=h/wune pulsation w 
non pas égale a la pseudo-pulsation w4, mais à celle pour laquelle le 
spectre des fréquences de l’oscillation amortie présente un pic marqué, 
à supposer qu'il existe. 

Ce spectre est, pour l’oscillation u(t) e-vt cos Wat, 


v + io 
w® + v2 — w? + 2ivo 


F(w) = Î e-({o+r)t cos (wat)dt = 
0 


d’où en module, 


(49) 


S(@) = |F(iw)| = | P+ ot ‘i 


(we + v? + w?)? — 4a?w? 
La courbe de S(w) présente bien un pic assez effilé lorsque € est 


suffisamment petit. Dans ce cas, on trouve dS/dw = 0 pour w = 0 et 
© = Om OÙ 


Om = N Ga aw? + 4v8 — v2 (50) 


Appliqué à notre cas, ce résultat fournit, avec ©? = k/m : 
Vv = c/2m = h2mom = nool2om (51) 
en méme temps que : 


Ga = Vw? — v2 (52) 
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Il suffit d’éliminer w » et v entre (50), (51) et (52) pour trouver wa 
en fonction de wo et de n. On trouve ainsi la relation : 


Maine) lo AGE (0 o°)] (53) 


ce qui, en posant wa = wo (1 — cin? — coy...) donne «1 = 1/8, cz = 1/32 
et donc une expression qui ne diverge de celle de Collar (44) qu’a partir 
du second ordre en »?2. On pouvait s’y attendre puisque w» se confond 
avec «4 au premier ordre près. 


6. SYSTÈMES A MÉMOIRE « COURTE } 


Le système représenté par le modèle (26), repris sous sa forme 
essentielle : 


mx + F’(|x|).x/4|x| + g@) = PO (54) 


s’il a le mérite de la simplicité, a cependant le tort de ne pas donner un 
cycle d’hystérésis du type habituel (fig. 5a pour « — 0). 


Fig. 5 
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On remarque qu'il suffirait d’inverser le rôle des sommets de la 
courbe avec celui de l’origine pour passer de la courbe a à la courbe b 
(fig. 5), cette dernière étant semblable à celle de la fig. 1. 

Or ceci peut s’obtenir par l'introduction de ce que l’on pourrait 
appeler une mémoire courte, en ce sens que l'influence du terme non- 
linéaire de l'équation précédente se ferait sentir à partir du dernier 
extremum atteint par la variable x. Ceci demande évidemment un ajuste- 
ment des autres termes et l’on obtient ainsi : 
nx Spl) F’(|x| = |Xm ) Aix Gr 2(x) — 2(Xm) = P(t) — P(tm) (55) 
étant entendu que xm désigne le dernier extrémé atteint. En fait, ces 
extrémés s’obtiendront pour x = 0 et à ce moment la fonction géné- 
ralisée x/|X| devient aussi nulle (moyenne entre les valeurs — 1 et + 1 
a gauche et a droite). Contentons-nous ici de quelques indications 
élémentaires concernant le cas schématique d’une linéarité fragmentaire, 
pour lequel (55) devient, avec F(a) = 2k8a? 


x 


pi 


mux = Xm) = kB(Ix| = |Xml) = & K(x — Xm) = P(t) = P(tm) (56) 


Fig. 6 


ge 


Le probleme des oscillations naturelles amorties doit étre traité 
en posant P(t) = 0, soit alors ao l’amplitude initiale, avec x9 — 0, et 
Xmo = 0. Les équations se présentent alors de la façon suivante (fig. 6) : 
pour Of << fom avec 0er) 0, 


mxX + kB(\x| — ao) + k(x — ao) = — P(O) = — kao (57) 


d’où, en intégrant après multiplication par x, et en prenant x (4) = 0 
et C2) — do: 


mx? — k( +8) ee. — 2k Ba- 
(58) 
—mx2 + KA — B) a? = — 2kBaoas 
COLE 
x2 = OR (1 —6) (59) 


et par addition : 
a2(1 —B) — ai +f) = — 26ao (ao + a2) (60) 
Ensuite, pour fg < t < ta, avec x > 0 et X(44) — 0, en notant que 


cette fois le second membre de (56) est nul, mais qu’au premier membre 
— MX¥m = kB (az — ao) — kag d’après (57) : 


mx — kB (|x| a2) + k(x + a2) =k [a2 + 8 (ao — a2)] (61) 


d’où, avec x(t3) = 0 et x(t4) = aa, le bilan énergétique des intervalles 
correspondants : 


mx? — Kk(1 + B) à = 2k (ao — 2a2) ae 


(62) 
— mx + kl — 6) at = 2k (ao — 2a2) aa 
On en déduit la valeur de xg, et par addition : 
a? (1 —6) — ae (1 + 8) = 28 (ao — 2az2) (az + aa) (63) 
De (60) et (63) on déduit, lorsque @ < 1: 
aan 10/2 
(64 
aalaz = 1 — 6 + 362/2 
ce qui indique un taux de décroissance : 
aalao = 128-262. (65) 


comparable jusqu’au second ordre avec celui (47) observé dans le 
système sans mémoire « courte ». 
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REMARQUE 


Dans le cas de l’hystérésis magnétique, on peut songer à construire 
un modéle mathématique suivant les mémes principes, mais en notant 
qu’ici le terme d’inertie fait défaut (m = 0). Un modèle qui apparem- 
ment se rapproche des conditions observées expérimentalement (avec 
toutefois un écart sur la valeur des pertes cycliques données par Stein- 
metz), serait 


=i (ae = Bali) À + k1(B — Bm) + ks(B5 — B3) = H — Hm 
‘B 

H et B désignant respectivement le champ et l'induction: le rôle du 
terme € > 0 consisterait à effacer la discontinuité de pente qui, sans 
lui, apparaitrait aux points B = O0 du cycle. 

Nous reviendrons sur ce cas ainsi que sur le probléme plus général 
des oscillations entretenues, dans une analyse des résultats fournis par 
des systèmes analogiques. 


ANNEXES 


1. On sait que du point de vue des mathématiques appliquées 
le maniement de certaines fonctions telles que celle de Dirac, utile 
dans l’analyse harmonique, a été grandement facilité par les travaux de 
Temple [6] et de Lighthill [7]. 

Rappelons que Lighthill considère des fonctions f(x) partout 
continues ainsi que toutes leurs dérivées, dont la limite pour |x| infini 
est au plus de l’ordre de x, d’une part, pour un N fini donné s’il s’agit 
de la classe des fonctions « passables » (fairly good), d’autre part pour N 
quelconque (même négatif) s’il s'agit de la sous-classe des « bonnes » 
(good) fonctions. 

Les fonctions généralisées sont alors définies par une séquence 
infinie de bonnes fonctions fa(x) telle que, F(x) étant une fonction 
passable quelconque, la séquence FA JnFdx, converge uniformément 
vers une limite finie. a 

Les indices ne peuvent correspondre aux valeurs successives attri- 
buées à un paramètre, et par exemple la fonction de Dirac d(x) est une 
fonction généralisée qui peut être définie par la condition que cette 
limite soit 


+00 
(3 3(x) F(x) dx = F(0) (a) 


— 0 


pour toute fonction passable F(x). 
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En fait les fonctions généralisées définies de cette façon ne doivent 
jouer un rôle que sous des intégrales, où on peut leur appliquer les 
régles valables pour les fonctions continues, puisque tous les membres 
des séquences qui les définissent sont de « bonnes » fonctions. 

C’est ainsi que la dérivée d’(x) de Ô(x) peut se définir par toute 
séquence donnant la limite : 


+00 +0 
i d'(x) F(x) dx = 8(x) F(x) | — 


CO 


+00 

if d(x) f(x) dx = — F'(0)  (b) 
On remarquera aussi que la somme ou le produit d’une fonction généra- 
lisée et d’une bonne fonction donnée, reste une fonction généralisée. 

Dans notre cas, si l’on admet que 2% (x, x) appartient à la 
classe des fonctions généralisées en x, il est évident que les seules fonctions 
«singulières » qui doivent nous intéresser sont celles qui se manifestent 
en x = 0, c.-à-d. (x), 5/(X), etc., puisque la dissipation cyclique est 
indépendante de la fréquence w > 0, et doit apparaître déjà aux vitesses 
les plus faibles, pour n’importe quelle amplitude. 

Il ne reste donc plus qu’à écrire 2; sous la forme : 


DE = fox, x) + fix, X) S(O + f2(x, À) YO + fax, X) SX) +... © 


où les fo, fi etc. sont de bonnes fonctions en Xx. 
Reprenons alors (17) et (18) sous la forme : 


T/2 5 T/2 3 
f[ ga | 2", %dr=...=0 (d) 
x Le 
0 0 


(ces intégrales restant nulles pour une demi-période a cause de la symétrie 
du système), nous voyons que la série de Taylor, applicable à fo (x, x): 


2 
folx 0) = fold) — fé D+ fon O58) — (e) 


multipliée par x? et intégrée de 0 à T/2, respectivement dans les inter- 
valles où x > 0, et x < 0, en exceptant les points où * s’annule, donne : 


T/2 LED he 
i fo(x,0+) x dt = 0 — : fo(aë, 0+) ee dt 
0 0 


pour toute valeue de a, de sorte que fo (x, 0 +) = 0. 
De méme on a: 


Wy 


LANDE; Te, 
hae (x,0+) Se (x,x) Xi a Eu ign 02° 


15 


et donc : 
T/2 


1/2 > 
/ : ae. es 4 LA 3 
f'. (x, 0H) ¥8dt = o= |! fy. (aE, 0) E8 dr 
0 


0 


pour toute amplitude a, d’ou : 
fi. (x~,0+) =0 


x 


On trouve de méme : 
A ty (EX, ) = 0 


et ainsi de suite, de sorte que la « bonne » fonction fo(x, X) est nulle sur 
tout l'intervalle 0 à T/2. 

Enfin, la première condition (d) montre que les fonctions générali- 
sées 5’’(X), 3/’(X), etc. ne peuvent subsister dans @; que multipliées 
par des fonctions de la forme XF3(x, x), X?F4 (x, X), etc. 

En effet, on a : 


EL En e\ O77 72 25) 2 . 
2 falx,%) 8°"? dt = Î fa (x,%) 3” (8) 2/5 di 
1 len 


= fod" ue |” - Hi "SOU, 8S) +L, CX) + 28/8) 


x ‘ 
m x m 


— fi. (X2°X/X2)] dx 


= Î OU, GR + Le, D + 2 (IR) XIE] di 


=f; (4,0) + 2fs (a,0)/Xa, etc. mn 


ce qui demande /3(x,0) — 0 pour tout x, c.à.d. fg (x, *) = XFs3 (x, X) 
etc., sinon il apparaîtrait des dérivées supérieures Xa, Xa, etc. (dont la 
valeur dépend des forces extérieures et qui ne peuvent donc caractéri- 
ser des paramètres intrinsèques du système, tels que les facteurs de 
dissipation). 

Dans ces conditions l'intégration de @; ne fait pas apparaître, 
par les termes XF38’(x), X2F45’’’(X) etc, de contributions formellement 
nouvelles par rapport à celle fournie par les termes précédents, fi 3(X) 


et f25(X), termes auxquels nous pourrons donc limiter l’expression 
de g:. 


2. L'intégration de (21) se fait le plus commodément en partant de 
x — 0 et en superposant les intégrales valables respectivement pour 
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X > 0 et pour x < 0; de toute façon, l’addition d’une fonction d’inté- 
gration g(x) encore arbitraire a ce stade justifie ce procédé, qui, a notre 
point de vue, a l’avantage d'être symétrique par rapport à x. On a, 
en remarquant que dx/dx = x/X : 


x>0 r<0 
D (x, à) = co il ~u(—%) i ES +R GI IAE +e (0 
0 0— 


+ 


Di 


, «> x>0 
= u(x)[fi(x,0+)/2 + (fa(x,x) d(x) |] — fi SOLAS (XX) ol dx 


| 0+ 0+ 


T<0 <0 
— uW—X)[fi(x,0—)/2 + (foxx) 8(X)/] — if SOLS, CX) + fix 
10= —0 d * 
+ g(x). (g) 
ou la fonction échelon-unité u(x) est définie par 
MES) = ie d(x) dx 


— C0 


donc nulle pour x < 0, égale a 1/2 pour.x = 0, et à | pour x = 0. 
On obtient ainsi, en supposant f,, (x, O+) < © : 


B(x,x) = u(X)[C1/2)[AC,0+) —£ 40+] — 0+)8Q)] (h) 


ih — ENO Nia —) = F0—=)) — f2,0—) 00)] 2) 


Posons : 
fix) = hee (x,X) = 4h(x,x) (i) 


et notons que, avec : 
2h(x,0) = CO) A0) 
on peut écrire : 
2h(x,0+) = [A(X,0) + A(— x,0)] + 
Lx, 0) — A(— x,0)] + G0+) — h&,0—)] (3) 


et une formule analogue pour 2h (x, 0—). 
On peut maintenant transcrire (h) sous la forme : 


D (XX) = 2h(x,0)[u(x) — u(— X)] + [h(x,0 +) — A(x,0—)Iu(x) + u(— x)] 
— [F(x,0+) u(x) — F(x,0—) u (— HI 8%) + e@) (kK) 


ce qui, compte tenu de (29) nous donne bien l’expression (22) puisque : 


10) ue ey JA a) ue |. 
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Deuxième Section 


SCIENCES PHYSIQUES ET CHIMIQUES 


Sur les aberrations périaxiales de l'œil. III 


PAR 


A. BIOT 


SOMMAIRE 


On calcule les aberrations géométriques de l’œil théorique en la fovéa 
pour les radiations C et F. Puis, en tenant compte de la courbure des rayons, 
celles qui correspondent, pour la fovéa, à la loi de variation des indices dans 
le cristallin n? = n? + (3.1319325)y°, n, étant l'indice sur l’axe et y la distance 
à l’axe. On discute les résultats obtenus qui montrent en particulier que pour 
une valeur de la constante, dans la loi de variation de l’indice, égale à la 
moitié de celle qui est donnée ci-dessus on obtient une réduction au 1/10 des 
aberrations pour la raie D dans le plan focal, puis, pour C et F, une sensible 
superposition des «disques » de diffusion. 


1. Les aberrations géométriques de l'oeil théorique en la fovéa 
pour C et F. Nous avons considéré les mêmes rayons incidents paral- 
lèles que dans le travail I (1). Ils font un angle commun de 5° 00’ 347.3 
avec l’axe de l’œil et sont numérotés de | à 6. 

Nous avons donné, dans le travail cité, les indices adoptés pour la 
radiation F. Les indices pour C ont été calculés de la même manière 
que ceux pour F. Voici les valeurs adoptées 


cornée 1.37501 
hum. aq. 133553 
crist. 1.41741 
vitré 1.33414 


Le tableau (1) ci-dessous renferme les valeurs des hauteurs h; 
d’intersection dans le plan focal pour C, D, F, les valeurs pour D et F 
étant reprises du travail I pour la comparaison. Nous avons indiqué 
aussi les hauteurs A; d’intersection sur la cornée. 

Le rayon 3, calculé avec les indices relatifs à la raie D, passe par 
le centre de l'iris. Calculé avec les indices relatifs à C et F il n’en est 


() Ann. Soc. Scient. Brux., 74, 175, 1960. 
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plus de même. Toutefois les hauteurs h; relatives a des rayons qui 
passent, pour C et F, par le centre de l'iris coincident pratiquement 
avec les valeurs données ci-dessous. 


rayon Rs | 
© D | F 
1 + 1,75 1,4018 * 1,3897 1,3586 
2 + 0,80 1,4517 1,4517 1,4278 
3 — 0,2657 1,4625 1,4609 1,4568 
4 — 0,6815 14627" 1,4629 1,4634 
5 — 2,00 1,4885 1,4949 1,5108 
6 — 2,26 1,4995 1,5069 1,5267 


Remarquons encore que les valeurs Aj¢ marquées d’un astérisque 
n’ont pas été obtenues directement. Ayant en effet obtenu les autres 
nous avons remarqué que pour celles-ci le rapport |hic — hipl/lhir — 
hip! est égal à 0,395. nous avons admis qu'il en est de même pour les 
rayons | et 4. D'où les valeurs inscrites au tableau. 


2. Calcul de h; avec un cristallin à indice variable (C et F). Nous 
avons vu que, parmi les relations essayées, seule la formule n? = n°0 
+ by? fournit des résultats encourageants. Faut-il, pour C et F, adopter 
des valeurs de b différentes de celle que l’on a choisie pour D?? Nous 
n'en avons rien fait. Une variation (légère) de b introduirait une variation 
de n qui serait en effet du second ordre puisque l’introduction de b ne 
produit déjà qu'une faible variation de indice. Nous avons en consé- 
quence fait nos calculs en prenant b = {3,1319325} = bn. 

On trouve alors, pour la raie F, les relations suivantes pour les 
rayons courbes dans le cristallin 


rayon 
3 x + 38,21389 =  {1,9493946} Logio (y + 4/}? + 7,223877) 
4 x — 49,71136 = — {1,9489237} Logio (y + 4/v2 + 10,46031) 
5 x — 72,72190 = — {1,9470175} Logio (y + 4/72 + 23,48938) 
6 x —75,11225 = — {1.9466835} Logo (y + 4/2 + 25,76092) 


(n? = 2,0341605 + {3,1319325} y?) 
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Si, pour ces rayons (b = b1) on calcule les valeurs Ajp on constate 
que le rapport de |Aig — hip| en rayons courbes dans le cristallin à 
|hig — hip| en rayons rectilignes dans le cristallin est le même pour 
tous les rayons et vaut 1,030. En appliquant cette règle aux rayons 1 et 2 
courbes dans le cristallin on peut facilement calculer les A; correspondants 
à l’aide des nombres du tableau (1). Si l’on opère de la sorte il convient 
de déterminer le rapport pour des rayons de n° élevé; ce sont ceux pour 
lesquels les angles qui interviennent, et par suite la précision, sont les 
plus grands. 

Pour la radiation C nous n’avons déterminé que les équations, dans 
le cristallin, des rayons 5 et 6. On trouve 


I 


5 x — 71,97972 = — {1,9443589 } Logio (y + 4/72 + 22,94058) 


6 x — 74,33810 = — {1,9440376} Logio (y + +/y? + 25,09929) 


(n? = 2,009051 + {3.1319325 }y) 


Nous avons remarqué en effet que pour ces deux rayons le rapport 
| hip — hic|/| Aig — hip | vaut aussi 0.395. D'où les valeurs hic (6 =b1) 
pour les rayons (c) 1, 2, 3 et 4. 

Le calcul direct des h; pour b = b; est plus précis que le calcul 
direct de ces mêmes quantités pour b = be = b,/2 car les petits nombres 
qui interviennent sont plus grands dans le premier cas que dans le second. 
C’est pour cela que nous avons effectué les calculs précédents pour 
DD 

Nous savons pourtant que la dispersion des valeurs de /; est plus 
petite pour b = bz que pour b = b; (n = np). Nous avons donc ramené 
les valeurs h; trouvées ci-dessus pour b = b1 aux valeurs correspondant 
a b = be par une remarque faite antérieurement et avons établi de la 
sorte le tableau (2) ci-dessous 


Valeurs de h; pour b = be 


rayon (E Diff. D Diff. F OnD 


1 149075 0,0129 | 1,4248 0,0312 | 1,3936 | + 19 10° 03°'.6 
22 1,4762* 0,0068 | 1,4694 O01 0. Wi 4520 le 22251078 
3 1,4661 * 0,0017 | 1,4644 0,0041 | 1,4603 | — 4°05’ 04'".9 
4 1,4572* | —0,0003 | 1,4575 | —0,0006 | 1,4581 | + 5°07’ 45"'.6 
5 1,4564 —0,0066 | 1,4626 | —0,0164 | 1,4795 | + 8°33’ 44.3 
6 1,4647 —0,0077 | 1,4724 | —0,0202 | 1,4925 | + 9°04’ 29°”.9 


a ——]_—_—_—_—_— 
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La derniére colonne renferme, pour la radiation D, la valeur de 
l'angle que font, dans le vitré, les différents rayons avec l’axe optique. 
Le signe plus devant w, indique que le rayon converge vers l’axe apres 
sa sortie du cristallin; le signe — qu'il s’en écarte. 


4. Discussion des valeurs de hy pour b = bg. On reportera avec 
avantage les valeurs de h; sur des verticales C, D, F, à une échelle con- 
venable et à partir d’une horizontale commune. D'un côté de ces droites 
on placera les valeurs pour b = 0, de l’autre les valeurs pour b = bp. 

a) On remarque ainsi facilement que les points relatifs à b = be 
sont bien plus serrés (mieux groupés) que pour b—0, le rayon | montrant 
cependant chaque fois une notable sous-correction. Si l’on ne tient pas 
compte de ce rayon on constate que pour D la dispersion des b; est 
environ 10 fois plus petite que pour b = 0. Il y a un gain notable pour 
C et pour F également. 

b) Ilest permis de se demander pourquoi l'augmentation de l'indice 
avec la distance à l’axe dans le cristallin diminue l’aberration alors 
qu'à première vue celle-ci devrait être augmentée en raison d’une 
possible augmentation de puissance marginale. La diminution de l’aber- 
ration provient de la courbure du rayon. Celui-ci possède en effet, 
puisque b est > 0, une courbure tournée vers les y croissants (y est > 0 
vers le haut dans la moitié supérieure du cristallin et vers le bas dans 
sa moitié inférieure). De la sorte, la tangente au rayon lumineux courbe 
au point de sortie du cristallin est moins inclinée sur l’axe que celle de 
l’entrée. Elle coupe l’axe optique plus loin que ne le ferait un rayon 
resté rectiligne (indice constant). Et cela correspond à une « correction ». 

c) Si l’on néglige toujours l'effet très partiel du rayon | et,éventuelle- 
ment celui du rayon 6 on constate que les disques de diffusion corres- 
pondant à C et à F sont sensiblement coaxés et renferment la plupart 
des rayons de l'ouverture. Cela est nécessaire pour que ne soient pas 
perçu l'effet des aberrations chromatiques longitudinale et latérale 
qui ne sont pas perçues dans la vision normale dans laquelle est « neu- 
tralisée » la lumière colorée « débordante ». 

d) On peut rechercher s'il existe un rayon incident qui, quel que 
soit l’indice, rencontre le plan focal à la même distance de l’axe. On 
voit de suite qu'il en existe pratiquement un, compris entre les rayons 
3 et 4, plus voisin de ce dernier que de l’autre. La hauteur d’intersection 
peut se déterminer à l’aide des résultats du tableau (2). Elle correspond 
à des différences nulles entre hic, hip et hig On trouve ainsi, par la 
méthode graphique (II, $5,a) que pour la valeur heo = — 0,6228 
la valeur commune de hj = 1,4584. Cette valeur ne correspond avec le 
centre de gravité d’aucun des « disques » de diffusion C, D ou F. On 
te qu'elle est assez différente de celle qui correspond au rayon 
central. 
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e) On remarquera sur la figure proposée que le diamétre du disque 
de diffusion change peu quand on passe de D à C (x = 656,3 mm). 
Nos conclusions relatives aux h;p sont donc valables près de D vers €, 
c'est-à-dire dans l’orangé. On sait que l’on admet que l’œil met au 
point sur une image orangée. Il est d’ailleurs possible, en vertu d’une 
remarque faite plus haut (§ 2) de trouver toutes les valeurs h; corres- 
pondant à une longueur d’onde comprise dans l'intervalle Ac, Ap ou 
un peu en dehors de cet intervalle à l’aide du calcul direct d’un seul 
rayon lumineux — après que l’on a bien entendu déterminé à l’aide 
d’une formule d’interpolation convenable la valeur des indices de 
réfraction pour la longueur d’onde choisie. 

On peut opérer graphiquement à partir du résultat obtenu par le 
calcul direct qui permet de déterminer le rapport | hip — hi, | /| hin— 
hig |. Pour cela on établit à une échelle convenable (p. ex. 1 u en 4 = 
1 mm sur les échelles) le schéma de verticales dont il est question au $ 4, 
la distance de ces verticales pour D et F étant par exemple de 100 mm. 
On garnit ces deux verticales des valeurs h; du tableau (2) et on réunit 
par des droites les points qui correspondent au même rayon. On mène 
alors une nouvelle verticale à une distance de la droite I donnée par le 
rapport trouvé ci-dessus. Les intersections de la nouvelle verticale avec 
les autres droites donne les valeurs des h;;. 


5. Conclusions générales. a) Si l’on suppose que l’indice, dans le 
cristallin. varie avec la distance à l’axe il est nécessaire de tenir compte, 
pour calculer les aberrations, de la courbure concomitante des rayons. 

b) La loi de variation n? — n? + by? dans laquelle b est voisin de 
{4.8309025} réduit à 1/10 de sa valeur le diamètre du disque de diffusion 
calculé dans les mêmes conditions pour l’œil théorique normal. 

c) Effectués avec la même valeur de b, les calculs pour les radiations 
C et F donnent dans le plan focal de l’œil théorique modifié des disques 
de diffusion sensiblement centrés. 

d) La meilleure valeur de b correspond à une variation faible 
de l’indice au bord de la partie optique du cristallin. On conclut à une 
grande sensibilité des observations aux défauts locaux du cristallin. 

e) La valeur d’une «correction » géométrique des aberrations 
ne peut se déterminer que si l’on utilise les résultats obtenus pour cal- 
culer le chemin optique aberrant au « point » image. C’est ce que nous 
nous proposons de faire dans un prochain travail. 
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Sur les aberrations périaxiales de l'œil. IV 
PAR 


A. BIOT 


SOMMAIRE 


En se basant sur des résultats précédemment obtenus on calcule les 
différentes aberrations géométriques de l’œil théorique de Y. Le Grand au 
point d’intersection avec le plan focal du rayon principal du faisceau qui 
forme l’image et on en déduit les chemins optiques aberrants correspondants. 
Supposant ensuite que l'indice dans le cristallin varie suivant la relation 


nz = 1,422 + {4,8309025 } y? 


dans laquelle y est la distance à l’axe on calcule l’aberration sphérique longitu- 
dinale et la coma ainsi que les chemins optiques aberrants correspondants qui 
sont bien plus petits que dans le premier cas. 


l. Oeil théorique normal. L’aberration sphérique. a) Le calcul 
donne pour des rayons incidents parallèles à l’axe optique et distants 
de cet axe de 1,2 et 2 mm les résultats suivants 


b—0 
hcorm. | h crist. ab. sph.A x’ | wn en deg. @n en rad. 
mm | mm 
— ae — a : 
152 1,06099 — 0,1955 3° 07’ 08’ 
2 1,76868 — 0,5554 50 15 22€ {2,99660 } 


On OU, plus loin, « est l’angle d’intersection avec l’axe du rayon dans 
le vitré. 

b) L’angle «w restant petit et le système optique normal, on obtiendra 
une représentation suffisante de l’aberration sphérique longitudinale 
Ax’ en posant 

Ax’ = aw?, (2) 
la constante a étant déterminée par les valeurs de Ax’ et de w au bord. 
On trouve 


a = — {1,81141} 
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Dans les mémes conditions le chemin optique aberrant nA, n étant 
Vindice du vitré égal à 1,336, vaut 
nA = — naw’4/4 
soit 
nA = + 1,336 x {1,81141} x { 5,98640 3/4 = + {3,32156} 
_  Exprimons cette longueur en longueurs d’onde Ap = 589,3 mu = 
{4,77034} mm. On a 
TiO 
avec « = + 3,56. 
Dans le plan de meilleure mise au point, situé à Ax’M/2 c’est-à-dire 
à 0,278 mm en avant du foyer, l’aberration est réduite au 1/4 de la 
valeur absolue de celle que nous venons d'indiquer. On a alors une 
aberration zonale telle que 
4 = — 0,89 
nA étant nul au bord. 
On est loin de la condition de Rayleygh qui s’imprime par 
|a| < 0,25 
c) On peut, avec la méme approximation que précédemment, 
calculer l’aberration Ax’ en fonction de la hauteur d’intersection sur 
l'iris ou, ce qui revient au même, sur le cristallin. On écrira alors 


Ax’ = ah? 
La valeur de a’ se détermine à l’aide des résultats précédents par 
— 05554 =a & 176808 


On trouve 
a’ = — {1, 24931}. 
D'où 
Ax’ = — {1,24931} #2 (3) 
Cette relation permet de calculer l’aberration sphérique longitu- 
dinale pour des rayons qui coupent le cristallin en des points dont la 
distance à l’axe est la même que pour les rayons | à 6 des calculs précé- 
dents. Ces distances sont reprises dans la ligne 1 du tableau (4). Les 
aberrations Ax’ correspondantes se trouvent sur la ligne 2 du méme 
tableau. L’inclinaison & de ces rayons avec l’axe optique peut se 
tirer, avec une approximation suffisante des nombres de la colonne 7 
du tableau (2) du travail III qui précède. Ces inclinaisons sont sur la 
ligne 3 du tableau ci-dessous. En multipliant une valeur quelconque 
de Ax’ par la tangente de l’angle w correspondant on obtient la valeur 
de l’aberration transversale (ligne 4). Cette aberration, convenablement 
combinée avec les A; donne les valeurs hj (ligne 5) corrigées de l’aber- 
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ration sphérique et dues seulement a l’astigmatisme, la courbure de 
champ et la coma. Toutes les longueurs sont exprimées en mm. 


2. Oeil théorique normal. Astigmatisme et courbure de champ en 
la fovéa. a) Les aberrations se limitant au 3e ordre, il est intéressant de 
calculer la courbure de Petzval de l’œil théorique. La formule à appli- 
quer est 


dans laquelle n, est l'indice de réfraction du vitré, r le rayon d’une des 
surfaces de l’œil, n et n’ les indices correspondants. On trouve 


Rp = — 17,232 mm 


Un calcul facile montre qu’à la hauteur de la fovéa cette surface se 
trouve à — 0,0610 mm (en arrière) du plan focal. 

b) Astigmatisme. Le calcul exact donne pour les distances As 
et Ar qui séparent, sur le rayon central, le plan focal des focales sagit- 
tales et tangentielles les valeurs 


As = — 0,0784 et Ar = — 0,1074 


L’astigmatisme vaut donc 
As — At = 0,290 mm 
valeur en accord avec celle qui a été donnée dans le travail I et avec celle 
de la distance du plan focal a la surface de Petzval sur le méme rayon 
central. 

c) Au degré d’approximation que nous considérons les rayons 
incidents parallèles 1 à 6 étudiés plus haut et situés dans un plan méri- 
dien devraient se couper sur la focale tangentielle si l’aberration sphérique 
n'existait pas. Or celle-ci a été corrigée dans les résultats de la dernière 
ligne du tableau (4). Nous pouvons corriger ceux-ci de l'effet de l’astig- 
matisme et de la courbure de champ comme nous l’avons fait pour 
l’aberration sphérique. Ces aberrations introduisent en effet une aber- 
ration transversale, dans le plan focal, égale à At X tgw, l'angle w 
étant, pour chaque hauteur d’intersection sur l'iris, égal à la valeur que 
l’on trouve dans la ligne 3 du tableau (4). Les valeurs h’;, corrigées de 
At x tgw, donnent des valeurs h’’; que l’on trouve dans la ligne 3 du 
tableau (5) et qui correspondent au seul effet de la coma. 

d) Le chemin optique aberrant dû à l’astigmatisme et évalué au 
point de percée du rayon central avec le plan focal vaut (1) 


2 
no?h; 
4 


nÀ = — (C — A cos 29) 


, ( A. MARÉCHAL. Imagerie géométrique et aberrations. Editions Revue 
d'Optique, Paris, 1952. 
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avec, Ch 71. — lo OueL 
C = (At + As)/h? A= (As — At)/h? 
c’est-à-dire 
na)? 


4 


nA = — 


((Ar + As) — (As — Ar) cos 2¢] 
ou 
nA = + 1,336 { 3,93776 } [0,1858 +- 0,0290 cos 29] 
L’angle © est compté dans l’ouverture du système optique. 


Son maximum correspond au maximum de nA qui est donc 
nA = 1,336 x {3,93776} x 0,2148/4 = {4,78354} 
c'est-à-dire 
nA = + 1,03 Ap 

3. Oeil théorique normal. Coma. L'effet géométrique de la coma est 
représenté par les valeurs de h’% du tableau (5). On voit facilement que 
la courbe (hi, her) est symétrique par rapport à l’axe, ce qui correspond 
bien à un effet de coma pure. La valeur maximum de 3h’; (en valeur 
absolue par rapport au rayon central correspond à 36h’; = — 0,0095 mm 
environ. 

Par ailleurs on a, dans le plan méridien où A est maximum, pour 
la coma pure 


nA = nbh; sin w 
St 
db" = 3bh sin? 
D’ot 


nh =n Sh” sin «3 


Il s’ensuit que le chemin optique aberrant dû a la seule coma 
vaut ici en prenant w = 5° 20’ comme pour l’astigmatisme 
nA = — 1,336 x 0,0095 x {2.96888} = {4.59529} 
ou encore 
nA = — 0,67 ÀD 


Remarquons encore que, pour les hauteurs de rayons parallèles 
à l’axe égales à 1,2 et 2, le calcul exact donne 


heo nx” h/sin © — f h/sin wo — f — Ax’ 
Le — 0,1955 — 0,233 — 0,0378 
2,0 — 0,5554 — 0,6616 — 0,1062 (6) 


La dernière colonne de ce tableau indique l'effet de la coma qui 
augmente avec l’ouverture. 
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4. Si l’on assemble les chemins optiques aberrants au point d’inter- 
section du rayon central avec le plan focal on trouve, pour le chemin 
optique aberrant total 


(nA) = ZnA = (3,56 + 1,03 — 0,67) AD = 3,92AD 


Ce chemin correspond aux deux points du bord de l'iris situés 
dans le plan méridien. 

Si l’on suppose que par un choix convenable de la mise au point 
on puisse réduire l’aberration totale au 1/4 de la valeur ci-dessus on a 
encore 


(An) = + 0,98 Ad 


Cette valeur ne correspond pas aux observations qui donnent 
(nA) 7 < 0,25 Ap. De plus, elle correspond a une sous-correction au 
point de vue géométrique alors que l'expérience indique une sur-cor- 
rection au moins pour une partie de l’ouverture. 


5. Oeil théorique avec indice variable dans le cristallin. Aberration 
sphérique. Nous avons calculé l’aberration sphérique pour les hauteurs 
d’intersection 1.2 et 2 sur la pupille d’entrée (ou sur la cornée), la 
valeur de b adoptée étant be = {4.8309025 }. 

Les équations des rayons dans le cristallin sont 


x — 66,11087 = — {2,0987484} Logio (y + 4/7? + 4,166855) 


x — 94,39094 = — {2,0981737} Logio (y + Vy? + 12,02382) 


Les valeurs trouvées pour Ax’ et les wn» de sortie sont 


iF Ax On en deg. @n en rad. 
122 + 0,2663 3° 02’ 40” {2,12138} (7) 
2,0 — 0,1702 50 11’ 48” {2,95360 } 


L’aberration est repliée en sens inverse du sens habituel de correc- 
tion et elle s’annule dans l’ouverture considérée pour une certaine valeur 
de w. Dans la représentation de Ax’ en fonction de w il faut faire inter- 
venir un terme d’ordre supérieur. En se servant du tableau (7) on trouve 


hx = {2,19705 } De {4,34497 } 4 (8) 

Ax’ passe par 0 pour w = {2.92604 } qui correspond à une hauteur 
d’intersection sur le cristallin Aer = 1,6595 et une hauteur d’intersection 
sur la pupille d’entrée (ou encore sur la cornée) Heo = 1,8765. Le pre- 


mier nombre s’obtient en exprimant Ax’ en fonction de la hauteur d’inter- 
section sur Viris (ou le cristallin) ce qui conduit à 


Ax’ = {1,60220} fer — {1,16244} hg 
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On trouve le second nombre en multipliant le premier par le facteur 
1/0,8844 = heo/her qui est sensiblement constant pour toute l’ouverture 
de la pupille. 

Nous supposerons, dans ce qui va suivre, que l’on restreint l’ouver- 
ture de la pupille à & = 2 x 1,8765 = 3,753 mm. 

L’eil est alors corrigé géométriquement au bord où le chemin 
optique aberrant maximum vaut alors, au foyer 


1 1 
nA = na’ (3 -- À = na’ 6/12 


6 
a’ étant égal à { 4,34497}. On obtient 
nÀ = — 1,336 X {4,34497 } {7,55624 } {2,92052 } 
= — {4,95785} = — 1,54 À5 


On sait toutefois (2) que dans le cas où une aberration sphérique 
géométrique de la forme (8) est corrigée au bord il existe un point sur 
l’axe où le chemin optique aberrant ne vaut — avec un signe opposé — 
que le 1/16 environ du chemin optique aberrant que l’on trouve au 
bord de la pupille (4) avec le terme de degré supérieur. Il vient ici 
dans ce cas 

na’ (6 
nA = — — = — {404532} = — 0,19À 
16 6 t } ' 

Remarquons encore qu'en raison du mode de correction il y a 

surcorrection géométrique dans toute l'étendue de la pupille. 


6. Oeil théorique avec cristallin à indice variable. Coma. Observons 
d’abord qu'en raison de l'intervention des aberrations du 5° ordre 
et peut être d’autres aberrations du 3° ordre que les aberrations classiques, 
il est impossible d’opérer ici la dissection des aberrations géométriques 
que nous avons opérée pour l'œil théorique normal. II faut aborder le 
problème autrement. 

Le calcul exact donne les résultats ci-dessous 


heo Nae h/sinw — f  h/sino — f— Ax sin © 
1,2 + 0,2663 + 0,3045 + 0,0382 {2,72520} x 
2,0 — 0,1702 — 0,2077 — 0,0375 {2,95701 } Oo 


En tenant compte de l’aberration d’ordre supérieur, on peut écrire 


h/sinw — f = f(b sin?« + b’ sintw) 


OU 
h/sinw — f = f(bw? + b’a4) = fh(w) 


. di A.E. CoNRADY. Applied optics and optical design, EL., Dover, New- 
ork. 
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(Les coefficients b et b’ ne sont évidemment pas ceux dont nous 
nous sommes servi plus haut). On trouve 


b = {0,90887} b’ = — {3.05946} 
D’ou 
fO(@) = {0.90887} a2 — {3.05946} a4 (10) 

Cette expression est égale a zéro pour w = {2.9247}, c’est-a-dire 
pour une valeur trés voisine de celle qui fait Ax’ = 0. 

De ce résultat et des nombres de la colonne 4 du tableau (9) on 
peut déduire que la courbe qui représente la coma reste très voisine de 
celle qui représente l’aberration sphérique. La coma, ici reste bien plus 
faible que dans le cas de l’œil théorique normal. 

On pourrait d’ailleurs étudier la valeur duchemin optique aber- 
rant en se servant des relations (8) et (10) introduites dans la formule 


nA coma = nhi sin wo [fo(w) — Ax1/f 


Nous n’avons pas estimé utile de le faire (§7). 


7. Oeil théorique avec cristallin à indice variable. Les autres aber- 
rations. Nous savons que nous devons tenir compte des aberrations 
du 5e ordre. On n’a pas, que nous sachions, développé la théorie de 
ces aberrations pour les milieux à indice variable. Mais nous pouvons 
supposer que les 9 aberrations classiques du Se ordre existent ici aussi. 
Il est dans ces conditions pratiquement impossible d’aller plus loin 
que nous l’avons fait. 

On peut toutefois raisonnablement espérer que les améliorations 
constatées sur l’aberration sphérique et la coma s’introduisent finale- 
ment dans l’imagerie. 


8. Conclusions. a) Nous avons montré que les aberrations de l’œil 
théorique normal sont, pour la radiation D, plus grande et d’une autre 
nature que celles qui sont observées sur l'œil réel, 

b) Si l’on suppose que l’indice, dans le cristallin, varie suivant la loi 


go = 142 {4,8309025 } be 


on obtient pour l’aberration sphérique et la coma des valeurs compa- 
tibles avec celles qui sont observées et de même nature que celles-ci. 
c) Il est remarquable qu'une très faible variation de l’indice rende 


compte des observations. 
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Principe d’une méthode d’alignement 


PAR 


A. BIOT 


SOMMAIRE 


On décrit une nouvelle méthode d’alignement qui utilise deux lunettes 
et un miroir plan. 


1. Le problème. Etant donnés deux points A et B qui peuvent être 
distants de 50 m et plus et un point quelconque € compris entre A et B, 
il s’agit de placer, avec une précision donnée, le point C sur la droite AB. 


2. La méthode. On dispose en A et B respectivement deux lunettes 
opposées, pointées l’une vers l’autre. Ces lunettes sont munies d’un 
réticule formé par exemple d’une croix de traits doubles. Immédiate- 
ment devant chaque lunette est placé un repère qui coïncide avec le 
point nodal extérieur de l'objectif correspondant. Cette dernière con- 
dition peut être réalisée mécaniquement. Ces repères sont amenés à 
coïncider avec les points A et B. 

On dispose, par ailleurs, d’une lame à faces bien planes et bien 
parallèles rendues toutes deux réfléchissantes. Nous appellerons cette 
lame le « miroir ». 

Pointons la lunette A sur le point B. Plaçons le miroir sensiblement 
au milieu de AB, ses faces étant à peu près perpendiculaires à AB. 
La lunette A observe maintenant par réflexion dans le miroir et, par 
réglage du miroir, peut être amenée à viser le repère A. Le miroir est 
alors perpendiculaire 4 AB. Réglons maintenant la lunette B par auto- 
collimation sur le miroir : sa ligne de visée est parallèle à AB. Les deux 
lunettes sont réglées et leur réglage ne devra plus étre modifié. 

Soit maintenant un point € quelconque situé entre A et B. Il s’agit 
de le placer sur la droite AB. Mettons le miroir entre C et B, sensiblement 
au milieu de l’intervalle CB. A l’aide de la lunette B disposons-le per- 
pendiculairement à AB. Pointons le repère placé en C à l’aide de la 
lunette A et déplaçons le repère jusqu’à ce que son image se fasse sur 
le réticule de la lunette A : le point C se trouve alors sur AB. 
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Le point C une fois en place, opérons avec les autres points D, etc., 
a régler, comme avec C. 


3. Exemple. Supposons AB = 50m. Prenons deux lunettes identi- 
ques dont l’objectif a un diamètre de 40 mmet dont la longueur focale 
vaut 400 mm. La précision de pointé d’une telle lunette — en air homo- 
gène — est d’environ 1/10 de sa résolution c’est-à-dire (3/10)’’ = 0,3’’. 
A 50m cet angle correspond à 0,3 x 50.000/200.000 mm = 0,075 mm. 


4. Remarques a) Les deux lunettes étant supposées identiques on 
peut, après avoir opéré comme il est dit plus haut avec les lunettes 
A et B, recommencer l’opération pour tous les points C, D, etc., en 
pointant la lunette B sur le point A et en réglant la lunette A par auto- 
collimation sur le miroir. 

b) Le non-parallélisme des faces du miroir peut introduire une 
erreur de pointé. Pour éviter cette erreur on doit faire attention à con- 
server à l’arête de l’angle des deux faces la même direction dans le plan 
perpendiculaire à AB lors des déplacements du miroir. On peut aussi 
faire sur la lame-miroir une surface bien plane qui sera le miroir et 
corriger l’autre face de manière à réaliser une lame a faces « optiquement 
parallèles ». L'erreur de pointage à craindre sera plus petite mais il 
sera toujours utile de prendre la précaution indiquée ci-dessus. 

c) Les caractéristiques optiques des lunettes à utiliser dépendent, 
dans une certaine mesure, de la distance AB. Toutefois, dans un genre 
d'applications donné, les distances AB peuvent en général être consi- 
dérées comme étant du même ordre. 

d) Dans tout cas particulier d’application la précision du pointage 
en un point quelconque de AB est la même qu'en B. 
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Modèles Statistiques dans l’étude de la réaction 
A127 (d,x) Mg? 
PAR 
G. DEMORTIER (*), J.C. JODOGNE et P.C. MACQ 


Centre de Physique Nucléaire — Université de Louvain — Belgique 


SOMMAIRE On détermine les paramètres statistiques dans l'étude de la 
réaction 


Al?? (dja) Mg? 
Dans les notations de T. Ericson (1), on trouve 
T =1:7 Mey «t=2 Mev 1 =18Mer a=5.3 Mey™ 


INTRODUCTION. La physique des états excités du noyau connait 
plusieurs voies d’approche. Depuis ces dernières années, la méthode 
statistique connaît un regain d'intérêt (1) grâce au développement de 
la spectroscopie tant des neutrons que des particules chargées; on a 
étudié les niveaux des noyaux par excitation et désexcitation jusqu’à 
10 Mev au-dessus du niveau fondamental. Partant de la notion d’en- 
sembles nucléaires excités en équilibre (le noyau composé), on se pro- 
pose d’en étudier les modes de désexcitation avec un minimum 
d’hypothèses sur les composantes individuelles de cet ensemble (les 
nucléons). 

Dans cette optique, nous avons étudié la désexcitation par parti- 
cules a du Si ?° excité par capture de deutérons dans Al??. Nous évitons 
ainsi toute perturbation due a des effets de couches. 


MÉTHODE EXPÉRIMENTALE. Le faisceau analysé de 11.85 Mev du 
cyclotron du Centre de Physique Nucléaire de Louvain est focalisé 


(*) Boursier L.LS.N. 
(1) T. ERICSON — Advances in Physics — 36 - 425 - 1960. 
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sur une cible de 270ugr d’Al (fig. 1). Les spectres d’énergie des part - 
cules « résultantes est analysé au moyen d’un détecteur, semi-conducteur 
de type p-n RCA Victor. Celui-ci est protégé des ions et électrons émis 
par la source par l’adjonction d’un champ magnétique placé entre la 
cible et le détecteur. Les signaux après amplification attaquent un 
analyseur RIDL à 400 canaux (fig. 2). 


vanne à rotule 


x —! 4 
| chambre de |L- = cible = 
= , 
Paden 5 = : == faisceau à 
| | 
RNA k | L «| 
J} 
= 
| 
—— 
echelle porte - ; 
absorbeurs += EE 


detecteur 


DS 3 a 
polarisation 


alimentation preampli 


Fig. 1 — Géométrie de l’expérience 


Les fentes d’entrée de l’aimant analyseur du faisceau incident 
avaient 1,5cm d'ouverture afin de rendre optimum la résolution du 
spectre a. Celle-ci était limitée à environ | % pour une énergie Ex de 
5 Mev. On a testé la linéarité du détecteur et calibré l’échelle d'énergie 
au moyen d’émetteurs «, Po 210 et Pb?1?, ainsi que par les réactions 


O18 + d_> NU + « 
C2 + d+B + « 


La limite inférieure du spectre utile est déterminée par le bruit de 
protons et deutérons; elle correspond au canal 105 de la figure 3 don- 


nant le spectre « relevé. 
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Fig. 2 — Schéma de l'électronique 


| 


Fig. 3 — Spectre 1 de la réaction Al?’ (d, x) Mg?5 


L’abscisse de cette figure est désignée par le numéro d’ordre des 
canaux. La correspondance en énergie est la suivante 


canal 200 + E, = 7.230 Mev 


Ex 
canal 


— 33,05) kev 


ANALYSE DES RÉSULTATS. La capture par Al?’ de deutérons d’une 
dizaine de Mey produit un noyau Si 29 excité à environ 29 Mev. La con- 
tribution d’interaction directe est très faible dans ce cas (2); cela peut 
s’expliquer par la hauteur de la barrière coulombienne (6.2 Mev pour 
ro = 1.5 Fermi (3)) : à part les premiers niveaux du Mg?® où l’inter- 
action directe est importante. il n’y a aucune inhibition coulombienne 
qui pourrait favoriser anormalement ce processus par rapport à la 
désexcitation d’un noyau composé. 

Partant de l'hypothèse où protons et neutrons forment un gaz 
de fermions libres contenu dans le volume V d’un noyau et possédant 
une énergie élevée (de l’ordre de 30 Mev), on montre (1) que le spectre 
des particules émises est décrit par : 


N(Eo) dEx © Ge (Ex) Exe (E) dEx 


exp {2+/aE} 
E2 
avec N(E,) : le nombre de particules « émises dans un intervalle 
dE, d'énergie autour de E,. 
o(E,) : la section efficace de formation du noyau composé a 
partir de particules « (3). 
E, : énergie de la voie «. 
E : l'énergie d’excitation du noyau résiduel. 
a : une constante proportionnelle à A de l’ordre de A/10 (4). 
o(E) : densité de uiveaux du noyau résiduel au voisinage de 
l’énergie E. 
On a porté à la figure 4 


R Ge (Ex) dEx 


N (Es) E? 
n —S ES 
G(Ex) Ey 


(2) N. CINDRO, M. CÆERINEO, A. STRZALKOWSKI — Nucl. Phys. 24 - 
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(5) BLATT et WEISSKOPF — Theoretical Nuclear Physics p. 352. 


(4) R.C. ALBERT, J. ANDERSON, C. WonG; P.R. 120 - 2149 - 1960; T. 
ERICSON (Réf. 1). 


il 


en fonction de E!’2; on en tire pour a (*) la valeur 


H =" 3.3 Mev 


+ + 
rh pa 


2 
g ra ene N(Ea)E | 
si Ex Te (Ex) 


+++ + __NEa ) 
Eux Oe (Eu) 


25 3 EM Mev?) 35 
Fig. 4 — Détermination graphique de a 


Partant de l’approximation d’un ensemble en équilibre d’entropie 
S, on peut étudier le spectre comme dans le cas d’un phénomène d’éva- 
poration résultant d’un apport d'énergie. On montre alors que 


N(Ex) dE RY Oc (Ex) Ex exp ( 3 Ee/T) dE, 


ou T est la température nucléaire du noyau résiduel définie par 


L . ds 
Ts dH, / Ex = Ey max 


(*) Notons qu’on définit souvent la distribution par 


N(E) dEx ~ Ge (Ex) exp {24/aE} dEx 
A partir de cette définition on aurait (fig. 4). 


a = 3.2 Mev-1 
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On a porté a la figure 5, 


= N(Eq) 
Ge(Ex) Ex 
en fonction de E,; on trouve 
T = 1.7 Mev 
100 


N(Ea) 
Fx Oe (Ex) 


10 


6 7 8 9 10 E {Mev) 


Fig. 5 — Détermination graphique de T 


Enfin définissons la température par 
NE) = te(E) 


où N(E) est le nombre de niveaux excités du noyau résiduel dans l’inter- 
valle d'énergie 0 — E. Si + dépend peu de l'énergie E, 7 est égal à T, 
la température nucléaire (°). 


(5) T. ErICsON — Nucl. Phys. 11 - 481 - 1959 


oh) 


100 


10 


7 E( Mev) 10 


Fig. 6 — Détermination graphique de f 
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Les résultats des mesures d’énergie des niveaux du Mg? jusqu’à 


7.62 Mev (8) nous ont permis de porter à la figure 6, N(E) en fonction 
de E; on en tire 


T = 2 Mev 


Remarquons que les valeurs de + et T sont très voisines de leur 
valeur asymptotique £, 


E,max \1/2 
t= | —— = 1.8 Mev 


a 


déduite de l’hypothèse d’un gaz de fermions (1). 


CONCLUSIONS. Une étude systématique des niveaux à partir de ce 
triple point de vue semble s'imposer. En-dessous d’une énergie E, de 
6 Mev il y a un net excès de particules « de basses énergies. Faut- 
il Vattribuer à la méthode expérimentale, à une barrière coulombienne 
surestimée (“) ou à une autre voie compétitive favorisée par une struc- 
ture en particules « 


Si?9 > n + Si8 + Mg*4*-+a-+n 
Des expériences sont en cours pour étayer ces conclusions. 


Novus tenons à exprimer notre gratitude au Dr. G. DECONNINCK 
pour diverses discussions, ainsi qu’au personnel du cyclotron pour son 
aide constante. 

Nos remerciements s’adressent également a Monsieur le Professeur 
M. de HEMPTINNE et à l’Institut Interuniversitaire des Sciences Nuclé- 
aires qui nous ont permis de mener a bien ce travail. 


(*) S. Hinps et R. MIDDLETON — Proceedings of the International Con- 
ference on Nuclear Structure — Kingston — 1960. 

(7) C. FuLMER et B. COHEN — Phys. Rev. 112 - 1672 - 1959. 
C. FuLMER et C. Goodman — Phys. Rev. 117 - 1339 - 1960. 
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Diffusion élastique de neutrons de 14 Mev aux 
petits angles 


PAR 


G. DECONNINCK (*), A. MARTEGANI (*) 
J.P. MEULDERS et J. STOQUART 


ABSTRACT 


Small angle elastic scattering cross sections are measured with 14 Mev 
neutrons for 4 groups of elements (Fe, Zn, Cu), (Cd, Sn), (Pb, Bi) and U***. 

Three effects may be responsible of an increase of the cross section in 
forward direction: magnetic coupling, polarisability and atomic electrons 
collisions. The results show an inexpected behaviour for Fe, Zn and Cu. 


La diffusion élastique des neutrons et des protons s’explique remar- 
quablement bien à l’aide du modèle optique de Bjorklund et Fernbach 
qui inclut dans le potentiel optique un terme de couplage spin orbite. 

Cependant, les mesures effectuées précédemment sur l’Aluminium 
aux petits angles (1) ainsi que celles d’Alexandrof sur U et Th (2) ont 
montré que l’écart par rapport aux prédictions du modèle peut devenir 
important lorsque l’angle de diffusion devient inférieur à 5 degrés. 

En fait, plusieurs contributions supplémentaires peuvent produire 
une intensification de la diffusion prévue par le modèle optique (3). 

1° L’interaction électromagnétique du neutron et du champ 


coulombien du noyau cible introduit un terme supplémentaire dans le 
potentiel 


eL 1 do 
——— — — 6. 
"2M2C2 r dr 
2° La polarisation du neutron par le champ coulombien introduit 
le terme 


Vpor = — 4a E® x étant la polarisabilité du neutron (a~ 10-42cm?2), 


Vem = — wp L 


(*) Chercheurs I.I.S.N, 
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3° L’interaction du neutron incident avec les électrons du cortège 
électronique peut produire une diffusion du neutron aux trés petits 
angles, a cause de la grande différence de masse entre les deux partenaires. 

Ces deux derniers effets sont encore mal connus et c’est la raison 
pour laquelle des mesures de diffusion élastique s’imposent aux petits 
angles. 

La polarisabilité du neutron dans le champ coulombien a été 
envisagée plus particulièrement par Alexandrov (2). 

Quoi qu'il en soit, ces effets contribuent à augmenter l’intensité 
de la diffusion élastique aux petits angles et l’effet sera d’autant plus 
important que la charge du noyau est élevée. On doit donc s’attendre 
à une augmentation plus importante dans le cas des noyaux lourds. 
Il ne faut pas perdre de vue que la contribution due au potentiel optique 
augmente avec Z également (suivant le facteur A?/3) de sorte que l’effet 
relatif sera moins important que l’on ne pourrait le croire à priori. 

Signalons enfin que suivant Barashenkov et Barbashov (3) l’effet 
dû à la polarisabilité du neutron diminue avec l’énergie. Les neutrons 
de 14 Mev utilisés dans nos expériences défavorisent donc la diffusion 
produite par la polarisation électrique du neutron. 


MESURES ET RÉSULTATS 


Grâce à la méthode décrite ailleurs (1) nous avons pu mesurer la 
diffusion élastique aux angles inférieurs à 5 degrés dans une série de 
corps groupés suivant les poids atomiques voisins 


Al (Fe Cu Zn) (Sn Cd) (Pb Bi U 


L’aluminium (voir ref. 1) montre un accord parfait avec le modèle 
optique sauf cependant en-dessous de 5 degrés cdm. Les résultats obtenus 
sur trois corps de poids atomique voisin de 55 (Fe CuZn, éléments natu- 
rels) montrent une curieuse tendance à la diminution de la section efficace 
lorsque l’angle diminue. Si l’on peut mettre en doute cet effet sur 
chaque corps pris séparément, il faut reconnaître que l’effet d’ensemble 
est manifeste. Certains points ont été repris à deux mois d'intervalle 
et ont donné des valeurs sensiblement les mêmes. 

Par contre, le Sn et Cd ainsi que Pb et Bi sont en accord avec la 


courbe théorique. 


(2) DECONNINCK, DEMORTIER et MARTEGANI. Ann. Soc. Sc. Brux. Tome 74, 
11 136, (1960): 

(?) ALEXANDROV JETP. 33 (1957) 294. 

(5) BARASHENKOV et BARBASHOV. Nucl. Phys. 9 (1959) 422. 


103 


Le cas de l'Uranium mériterait une attention particulière, car les 
seuls résultats probants sur la polarisabilité du neutron ont été obtenus 
par Alexandrof dans U et Th. 

Afin d’écarter la possibilité de détection de neutrons de fission 
le spectre diffusé a été observé à l’aide d’un analyseur à 400 canaux 
(RIDL) et les mesures ont été prises en utilisant une fenêtre (single 
channel) au lieu d’un discriminateur. 

La mesure obtenue est significative d’une augmentation de do/do el 
vers 3 degrés. 


CONCLUSIONS 


Une étude systématique de do/duw el aux petits angles révèle les 
phénomènes suivants : dans le cas de Fe, Cu, Zn une diminution vers 
l’avant, dans le cas de U une augmentation. 

Le comportement anormal de Fe, Cu, Zn peut difficilement s’inter- 
préter dans le cadre des théories exposées au début du présent article. 
Toutefois, un effet d’interférence entre l'amplitude de diffusion poten- 
tielle et les autres contributions, n’est pas à écarter. La méthode expéri- 
mentale elle-même pourrait être incriminée à cause de la détection de 
neutrons diffusés inélastiquement; cependant les sections efficaces 
de diffusion inélastique sont généralement d’un autre ordre de grandeur 
que les sections efficaces de diffusion élastique. 

Enfin il serait intéressant de poursuivre plus avant encore les mesures 
sur Uranium, afin d'observer une montée éventuelle de do/dw el 
comme celle observée par Alexandrov à 1 Mev. 

Nous remercions le Professeur M. de Hemptinne pour ses conseils 
et encouragements, ainsi que I’L.1.S.N. qui a permis ce travail. 
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